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Abstract 



In this article we give a detailed discussion of the mass perturbation theory 
of the massive Schwinger model. After discussing some general features and 
briefly reviewing the exact solution of the massless case, we compute the 
vacuum energy density of the massive model and some related quantities. We 
derive the Feynman rules of mass perturbation theory and discuss the exact 
ra-point functions with the help of the Dyson-Schwinger equations. Further 
we identify the stable and unstable bound states of the theory and compute 
some bound-state masses and decay widths. Finally we discuss scattering 
processes, where the resonances and particle production thresholds of the 
model are properly taken into account by our methods. 
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I. INTRODUCTION 

The massless and massive Schwinger models are two-dimensional QED with one mass- 
less or massive fermion. Both models have been subject to intensive study in the last 
decades. First of all, the exact solubility of the massless Schwinger model was discovered by 
J. Schwinger more than 30 years ago Later on, Lowenstein and Swieca constructed a 
complete operator solution of the massless model. By carefully investigating the role of large 
gauge transformations they found that these large gauge transformations have the effect of 
changing the vacuum ("instanton vacuum"), and therefore a superposition of all these vacua 
("9 vacuum") has to be introduced as a new, physical vacuum in order to render the theory 
sensible (requirement of vacuum clustering). In the usual terminology of gauge theories, the 
occurrence of instanton-like gauge-field configurations can be traced back to the fact that 
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the boundary of Euclidean space-time is a S 1 , and, further, that the first homotopy group 
of the gauge group is nontrivial, n 1 ([/(l)) = Z. 

In addition, the chiral anomaly and Schwinger terms are present in these models, and 
all these nontrivial field-theoretical features were one reason for the rising interest in the 
Schwinger model |jTT] 



31, mill- 



However, there is another reason for the study of the Schwinger model, namely its simi- 
larity to some aspects of QCD. All the above-mentioned features are present in QCD, too. 
In addition, a fermion condensate is formed both in QCD and in the massless and massive 
Schwinger models PJ3), |T5| , pT[ | , and confinement is realized in the latter models in a quite 
understandable way. In both models there are no fermions in the physical particle spec- 
trum. The lowest physical state is a massive mesonic fermion-antifermion bound state, the 
Schwinger boson. In the massless model this boson is free, and higher states are trivial free 
n-particle states. In the massive model these states turn into n-boson bound states, because 
the Schwinger boson is an interacting particle there |2l||2"g|J£l||Io| . 

When confinement properties are tested with external charges, the two models behave dif- 
ferently, too. In the massless model widely separated probe charges are completely screened 
via vacuum polarization, and the potential between the external charges approaches a con- 
stant. 

For the massive model this screening takes place as long as the probe charges g are 
integer multiples of the fundamental fermion charge e. When g ^ ne, the potential between 
the probe charges rises linearly for large distances. So " screening" is realized in the massless 
model, whereas true confinement takes place in the massive model (this feature was first 
discovered in |24| within a bosonization approach, and is further discussed e.g. in M,0,|39 



see also Section 6). Therefore, the Schwinger model was studied in order to get more insight 
into the phenomena of quark trapping and confinement ]TT|,|I5| , IE . 



Instanton physics was studied, too, in the Schwinger model |22| . |23|1 . 
In addition, because of its simplicity, the Schwinger model is often used for testing 

semi-classical methods 



some new methods of QFT, like light-cone quantization ||(| - |5T 



m , |33|j44| or lattice calculations |52[ - [|57 . 

As mentioned, the massless Schwinger model was first completely solved within the 
operator formalism ||, and the operator approach and the specific two-dimensional method 
of bosonization were mainly used in the subsequent years |T5| , p!6| , |2^j25| ,|9|| . It took some time 
until a path integral approach to the Schwinger model arose (mainly because of the problems 
with the nontrivial vacuum structure) JTIJ , [0 - , - , - |2C 



In this article all computations are based on the path integral formalism in flat Euclidean 
space-time. 

This article is organized as follows: 
After fixing our conventions we discuss some general features and the physical meaning of 
the vacuum angle 9. Then we turn to the massless model, because it will be the starting 
point for a mass perturbation theory. We discuss the relevance of the instanton vacuum 
and briefly review the exact path integral solution of the massless model to the accuracy we 
need in the sequel. 

In a next step we compute the vacuum functional and vacuum energy density of the 
massive model in mass perturbation theory by a method that is analogous to the cluster 
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expansion of statistical physics. Further we prove the IR finiteness of the mass perturbation 
theory and comment on the UV regularization. 

With the help of the vacuum energy density we compute the fermion and photon con- 
densates and the scalar and topological susceptibilities and generalize (to arbitrary order 
mass perturbation theory) some recent path integral computations |^4[ on screening and 
confinement in the massive Schwinger model. 

Then we develop the Feynman rules for the mass perturbation theory, which we need 
in the sequel. Actually, because of the chiral properties of the model these Feynman rules 
acquire a matrix structure. 

Further we compute the mass of the Schwinger boson by a direct application of this mass 
perturbation theory. 

From the equations of motion we derive, in Section 9, the Dyson- Schwinger equations 
of the model and use them to re-express the exact n-point functions in a way that is more 
suitable for the subsequent discussion of bound states, decays and scattering. From the 
exact two-point function we will be able to infer the complete bound-state spectrum of the 
model and find that it contains, in addition to the n-boson bound states, some further bound 
states (Section 10). 

Further, this exact two-point function provides information on the decay widths of all 
the unstable bound states. As an illustration, we will compute the masses and decay widths 
of the lowest bound states. 

At last, we discuss the scattering processes of the model, where we take properly into 
account all the resonances and higher production thresholds. All this may be derived from 
a resummed mass perturbation theory, without further approximations. 

We will use the following conventions, 

9iu> = V > e oi = -* > e nv£ vX = fli (!) 
7^75 = , ^Ihltilv = -2e Mi , , 75 = -r/V (2) 

H-»o) • '-(:;) ■ *-(;-0 < 3 > 

and we find it convenient to keep the notation (x , x\) in Euclidean space-time. As a 
consequence, the dual field strength pseudoscalar 

F= l -e, u F^ = -up (4) 

is imaginary (/? is defined in (6)). The Lagrangian of Euclidean QED 2 is 

L = $(i0-e4 + m)tf- ^F^F^ (5) 
and we use the following representation for the gauge field 

A^x) = -(dXz) + ^d u (3{x)). (6) 
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The pure gauge part, a, is unimportant for gauge invariant VEVs (as we deal with through- 
out the article) and may therefore be set equal to zero, a = (Lorentz condition). Using 
the electric charge e (which has the dimension of mass in two space-time dimensions) we 
define the Schwinger mass 

= 7 (7) 

71 

which is the mass of the Schwinger boson in the massless Schwinger model. Further we need 
the currents 

J M = # 7/i * , f; = ^75^ = W 
S = M , P = ^ 75 ^ 

S ± = ^P ± y , p ± = i(i± 75 ) (8) 
and define the Schwinger boson field, $, by 

J„ =: -^e^d"®. (9) 
V 71 " 

On the quantized theory the following equations of motion hold (more precisely, they hold 
on the physical subspace, see the remark after (103)), 

d^ v = eJ u (10) 
d^J b u — —F — 2imP, (11) 

7T 

which are the Maxwell and anomaly equations. In addition, the Dirac equation holds, but 
we do not need it. 

We use the following Feynman propagators: 

massless scalar propagator 

D (x) = ln(:r 2 - ie) , UD Q (x) = 5{x) (12) 

4"7T 



massless fermion propagator 



G °<*> = sfej =m (13) 



massive scalar propagator 
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D,(x) = -—KMx\) , £)» = __ (14) 



and Kq is the McDonald function with the properties 

Z 

Kq(z) — > —7 — In - + o(z 2 ) for z 



I 7r 

X (z) — > W — e~ 2 for z — > 00 (15) 

V 2^ 

where 7 is the Euler constant. 

In the sequel we will have to distinguish between two types of VEVs, namely VEVs 
with respect to the massless and massive Schwinger models. The former VEVs are always 
denoted by ()o, the latter ones sometimes by () m and sometimes without a subscript, (). 

Here we should add a last comment on our Euclidean conventions. They are chosen in 
such a way that they are as similar as possible to their 1+1 Minkowski-space counterparts. 
However, for this advantage we have to pay the prize that in the intermediate Euclidean 
computations some quantities are unphysical (e.g. imaginary field strength F and vacuum 
angle 9). 



II. THE 9 VACUUM 

Let us briefly recall the way the vacuum angle 9 enters the theory. The key observation 
is the existence of large gauge transformations Gi that do not leave the vacuum invariant 
and therefore give rise to an infinite number of vacua, 

Gi|0)=:|l) , GyO) = (G 1 ) n \0) =: \n) (16) 

and n is restricted to integer values by the requirement that gauge transformations must 
act uniquely on the matter fields (here the fermion). As a consequence, gauge fields may 
exist that tend to different pure gauges G n± for the time t — > ±00. They have instanton 
number k = n + — n_, and this instanton number may be computed from the Pontryagin 
index density u(x), 



v 



= J dxv(x) = J dxF(x) = k. (17) 



These instantons induce transitions between different, gauge equivalent vacua (16). There- 
fore, a superposition of these vacua has to be introduced as a new, physical vacuum, 



00 



\9)= £ e M \n) , G n \6) = e-»*\e) (18) 



n=— 00 



00 

(9'\6\9) = 5(9-9') e tke (0\6\k) (19) 

k=— 00 

and in (19) we used (m\0\n) = (0|O|n — to), which holds for gauge invariant operators O. 
Because all 9 vacua are invariant with respect to large gauge transformations, up to a phase, 
(18), the angle 9 has to be included as a new physical parameter into the theory. 
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From (19) we find that within the path integral approach the vacuum angle 9 may be 
included, e.g. for the vacuum functional of the massive Schwinger model, like 



oo „ 

Z(m,9)= jT e ike DA»D^D^e$ dx[ ^-^+ m )*-\ F ^ F ^ (20) 

k=— oo 

oo „ 

]T / DA^DtyDtyel^^-^™^ 1 ^ 2 ^^ (21) 



k=— oo ' 



(where we used (17)) and has instanton number k. Due to the anomaly there is a third 
representation for Z(m, 8). By performing a constant chiral transformation on the fermions, 

^^ e ^7s^r ) #^*j> e ^7 5j (22) 
the change of the fermionic path integral causes the anomaly, 

A = --P [ dxF, (23) 

7T J 

and in the action only the mass term is changed by this transformation. Choosing f3 — — |, 
the OF term is cancelled and we find 

Z(m,6) = DA^D^D^eJ dx ^^' e ^ +mcos8 ^^ +imsme ^ 5 ^ + ^ F ^ (24) 

k=— oo 

Rewriting the ^-dependent part of (24) like 

m(cos OS + i sin BP) = m{e ie S + + e~ i0 S-) (25) 

we conclude that general VEVs will not be holomorphic in me %e but depend on its complex 
conjugate, too, () m (me !9 ,me" !9 ). 



There is a nice physical interpretation of 9 that was extensively discussed in pq| , which 
we want to present now. Ignoring the fermions at the moment and treating F as the 
fundamental field in (21), we find the equation of motion 

F = --. (26) 

Z7T 

So 6 may be interpreted as a constant background electric field. 

[Remark. There is a simple and very instructive way of discussing this feature in the 
Schrodinger picture, which we want to describe briefly (this was shown to us by R. Jackiw). 
The Hamiltonian of the theory without fermion reads, in the gauge Aq = (remember that 
Aq has no conjugate momentum; F = OqAi] we temporarily introduce a finite length L for 
the space direction X\) 



1 r L 

H = - / d Xl F 2 

2 Jo 



The quantum theory is described in the Schrodinger picture by operators Ai(x) and F(x) = 
ij-^j^j acting on wave functionals ^[Ai]. The Gauss law for physical states simply reads 
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diFfrWAt] = dn— — = 

with the general solution 

#L4i] = /(a) , a := [ dxiA^x) 

Jo 

where / is an arbitrary function. So there remains only one degree of freedom (the zero 
Fourier component of A\). The Schrodinger equation reads (E . . .energy) 

1 r L 1 r L 5 2 

- / dx 1 F 2 (x)f(a) = -- / dxi 2 f(a) = Ef(a) 

2 Jo 2 Jo oAf{x) 

with the solution 

f( a ) = e~ iFa , E = / d Xl — = L — , 
JK J ' Jo 2 2 ' 

where F is the eigenvalue of F, 

Ff(a) = (F)f(a) =: Ff(a) 

So F is just a constant field with energy density ^y. Finally, when we perform a large gauge 
transformation eA\ — > eAi — d\\ X(L) = A(0) + 2nn, ^f[Ai] changes according to 

1 



- -5iA] = *[Ai]e _ 



Therefore, the state vector ^[-Ai] changes under large gauge transformations precisely like 
the 6 vacuum (18), provided we make the identification F = — equ. (26), which again 
shows that 9 is a constant background field.] 

When matter is included, it was shown in [^J that whenever 9 > n (9 < — tt), it is 
energetically favourable to create a real particle-antiparticle pair (with fundamental charge 
±e) that partially screens the background field (26). This explains the angular character of 
9. 

Quantum effects further change the picture. For the massless Schwinger model (m = 0) 
it is obvious from (24) that the background field is completely screened. The anomaly (23) 
stems from the vacuum polarization graph in QED2 (for a lengthy discussion and explicit 
computation of the QED2-anomaly see e.g. |I7 JT9[| ), therefore the screening is due to vacuum 
polarization. In the massive case this screening is not complete and we will find (see Section 
6) 

(F) m m sin 9 + o(m 2 ). (27) 
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III. INSTANTON VACUUM AND ZERO MODES 



The following two sections are devoted to a short review of some properties of the massless 
Schwinger model which we need in the sequel. First we remember that the exact fermion 
propagator in an external gauge field may be given explicitly, 

G p {x, y) = e^^-^Goix - y) (28) 

(more precisely, (28) is correct for gauge fields with vanishing instanton number, k = 0. 
For k G 13 acquires an additional term. This term, however, does not contribute to the 
path integral computations which we want to perform and may therefore be omitted, see 
the remark after (32)). 

Further we should recapitulate some wellknown properties of the Dirac operator in an 
instanton field. For sufficiently simple space-time manifolds the Dirac operator 

ty=$ + ie4 (29) 

in a gauge field with instanton number k has precisely \k\ zero modes, which have positive 
chirality for k > and negative chirality otherwise. For the Schwinger model these zero 
modes may be computed explicitly (see 



*f fc = -^L(x-)V A ( 1 ) , k>0 , i = O...Jfe-l 



'27T ' \0, 



K k = -^{x + ) i0 e- t(3k r j , k<0 , i„ = 0...|Jfe|-l, (30) 



'2ti 



x = X\ + ixo , x = X\ — ZXq. 

Next let us investigate the pure fermionic part of the path integral in an external gauge field 
(for a proper treatment of the zero modes we keep a small fermion mass m that we set equal 
to zero at the end of the computation). 
Because of the identity 



Z[A k u \ = lim / D*D*ef* B *W* n! ) 9 



= lim detUB + m) = lim m |fc| det'(i0+ m) (31) 

the vacuum functional obviously vanishes for k ^ (the prime indicates omission of the zero 
eigenvalues) . 

This remains true for VEVs of operators containing only gauge fields, because they 
do not influence the fermionic integration. On the other hand, fermionic VEVs may get 
contributions from nontrivial instanton sectors. 

By the use of Grassmann integration rules a general fermionic VEV may be written like 



(^{x 1 )^{y 1 ) . . . ^ a "(x n )^(y n )) [A 
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l im V (-) jW nE ^ ( 1^" )(l/T(i)) m^TJ'iiM + m). (32) 

Here ipi(x) and A; denote eigenfunctions and eigenvalues of the Dirac operator, and the sum 
just means summation of all possible Wick contractions (for a more detailed discussion see 
e.g. lull). 

This expression may give a nonvanishing result when in precisely \k\ Green functions 
only zero modes contribute, because then the m' fc ' factor is cancelled by a factor m~' fc '. Of 
course, every zero mode has to occur exactly once because of the Pauli principle (the Pauli 
principle is not explicitly written down in (32); however, because of the permutation sign 
factor (— ) CT ( 7r ) there is a pair-wise cancellation of all terms with identical eigenfunctions). 
Therefore, a VEV of n fermion bilinears, like (32), may get contributions from instanton 
sectors k = —n, . . . , n. 

[Remark, for a nontrivial instanton sector k ^ 0, the remaining Green functions in (32) 
(the exact fermion propagators (28)) in principle should be constructed on the subspace that 
is orthogonal to the zero modes, i.e. they should satisfy 

1*1-1 

^G^(x,y) = lS(x -y)-J2 (v) 

io=0 

and, consequently, read 

G Pk {x,y) = e^ 5Mx) G (x - y) e ^ A ^ - e ^ rMx) £ c ioX i k (x)^ k (y) 

io=0 

1*1-1 

= e^ 5 ^ {x) G (x - yy^ 5 ^ - J2 c i0 ^ k Xi k (y)e il5Mv) 

io=0 



where 



hxi (x) = e *wAW*g- (x) , d»xl k (vh, = K k (y)e* 



(the Cj are related to the zero mode normalizations and are unimportant for our argument). 
The important point is, of course, that the additional term for G@ k contains the zero modes 
and, therefore, does not contribute to the fermionic path integral (32) due to the Pauli 
principle. As a consequence, in our computations we may use the simple form (28) of the 
exact fermion propagator for all instanton sectors (this argument can be found e.g. in pip.] 
For further conclusions we have to specify the types of fermionic bilinears. E.g. pure 
vector-like VEVs get contributions only from the trivial sector k = for the following reason: 
the zero modes contribute like ^i ^i ~ P+ (for k > 0), the exact propagators are given in 
(28), and therefore all vector-like VEVs look like (for k > 0) 

trP +7M1 G^ 2 Gf 7^3 • • • = 
trP+Tw-P+Tw^sTwi • • • = tlP + p ~l^^ • • • = , etc. 
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and therefore vanish. 

On the other hand, VEVs of scalar or chiral bilinears do get contributions from the k ^ 
sectors. 

For densities, like S(x) = ^f(x)^(x) (which are the only objects we need in the sequel), 
this may be inferred in a very easy fashion from the various representations of the vac- 
uum functional, (20) - (24). Suppose for the moment that the fermion mass is space-time 
dependent. Then the scalar VEVs of the massless model are given by (see (20)) 

n 1 n S n 

(II ^)*(^)>o = — — I] - r --Z[m{x), 6]\ m=0 = (U(e l9 S + (^) + e^S.(xM = ° 
i=i ^(U, tf) i=l om{Xi) j=1 



= e^CQ ^ + (^))o =0 + e l( "" 2)e £<S-(*i) II S+fa))** + ■■■ + e- e <n (33) 

i=l j=l i =1 i=l 

where we used (24), (25) in the first line. 

Further we know from (20) that a factor e lkd indicates that the term stems from the 
instanton sector k. Therefore, we may draw the conclusion that for a n-scalar VEV the 
sectors k — n, n — 2, . . . , — n contribute. In addition, we find that for a VEV of n + positive 
chirality densities S + and n_ negative chirality densities S*_ only the sector k = n + — n_ 
contributes and that "instanton number equals chirality" . 

[Remark: we gave a quite explicit construction of the vacuum structure of the massless 
Schwinger model, because we need it for our further calculations. However, if one is only 
interested in the vacuum structure itself, it may be shown to be an almost trivial consequence 
of gauge invariance. First, imposing the Lorentz gauge condition a = in the representation 
(6) of the gauge field does not uniquely fix the gauge. There remains a residual gauge freedom 
to introduce functions (3 that fulfill the condition Uf3 = 0. Of course, a constant ft = c is a 
possible choice. Requiring gauge invariance and using the anomaly result (23) (for a constant 
chiral rotation c) we find for the vacuum functional 

J D^D^Dfte 3 ^^^ = J D^D^Dfte 8 ^ 5 ^^ = J D^D^D(3e s ^^~ 2tc ^ I ' dxF 

and conclude v — ^ / dxF = 0. Therefore, for the vacuum functional only the sector 
k = may contribute. The conclusion remains the same for gauge field VEVs and for vector 
current VEVs. For scalar and chiral VEVs, however, things change. E.g. for the positive 
chiral VEV (S + (x))q we find 

J D^D^Dft^{x)P + ^{x)e s ^^' f3] = J D^D^Dft^(x)e inh P + e ic ^^(x)e s[ *^^- 2icu 

and conclude (75-P+ = -P+75 = P+) that u = 1. Therefore, here only the sector k = 
1 may contribute. The generalization to higher VEVs is straight forward, and we may 
indeed conclude that the vacuum structure of the Schwinger model is a consequence of 
gauge invariance. (For the massive model the same gauge invariance requirement leads to 
the equivalence of the different representations (21), (24) for the vacuum functional Z(m, 9). 
A similar discussion can be found in [26], and more about instantons and 9 vacua e.g. in 



■)] 
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IV. SOLUTION OF THE MASSLESS MODEL 



The key observation for the exact solution of the massless model is the fact that the 
interaction term in the fermionic Lagrangian may be transformed away by a chiral rotation 

= - e4)^ = ^e-^Htye-^^. (34) 

In the fermionic path integral such a chiral rotation causes the chiral anomaly. In the model 
at hand this anomaly may be computed for finite chiral rotations, too, leading to the result 
for general VEVs (for 6 = 0; iV is the path integral normalization) 

N J D^D^D(30(^,^,(3)e s = 
N J D^ f D^ f D(3eI dx ^^f6(e l ^ 5 ^ f , ^/ f e 1 ^ 5 , (3)e ScS (35) 

where 

1 ^ 



D = -(□-/$ , G(x)=ir(D lio (x)-D (x)) , T> x G{x - y) = 5{x - y). (37) 

The first term in S e s is the "photon" kinetic energy, the second one stems from the anomaly; 
f is the free fermion spinor. 
For a further evaluation the presence of zero modes for k ^ has to be taken into account. 
In 0{^,^,(3) all Wick contractions among the fermions \P = e l/375 \l// have to be done and 
the corresponding number of k zero modes and propagators have to 

be inserted, according to our discussion in the last section. All this may be written down in 
a compact way by introducing the generating functional for fermions in the instanton sector 

k—l 

Z k [fj,n] = e ike N f Df3 k J] (fj\*%)(*%\ri)e- i f^W& k ( x > y ^e s '*, (38) 

«o=0 

where 77, fj are Grassmann-valued external spinor sources. 

Both exact propagator (28) and zero modes (30) depend exponentially on /3, therefore 
they add linear (3 terms to S c g, rendering thereby the (3 path integral Gaussian. As a 
consequence, all VEVs may be computed explicitly, as we now demonstrate briefly. 

E.g. for the chiral density (S + (x)) only the k — l sector contributes; we have to insert 
one zero mode (30) and find 

(S+(x)) = e w N [ DptiP+^jWWWp+ef* 1 ^? = 
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— e 2G (°) = e j9 ^ =: e ie - (39) 
2vr 4vr 2 v ; 

where we completed the square and performed the (3 integration in the first step; further we 
introduced the fermion condensate 

£ ee mxMx))t° = (40) 

Analogously we find for the two-point functions 

(S+(x)S-(y)) — —N [ D(3eI dz W- 



•tie^W-PWp+Goix - y)P~G (y - x) 



p 4(G(0)-G(*-v)) 7 

4tt 2 (x - y) 2 Mvr J 1 J 



(here only one of the two possible Wick contractions contributed due to trP + G = 0) and 

; ~ V) 2 c 4(G(0)+G(x-y)) = c 2ig/ e V 
47T 2 47T 



(5 + (x)S + (i/))o = . . . = e ««&j£ e ^(o)+G(*-i()) = ^(fZ^jV^**-*). (42) 



(the details of all these computations can be found e.g. in |18| , [19| ). 

Observe that in both cases the massless propagator part of G(x — y) is cancelled by a 
contribution steming from the free fermion propagator or from the zero modes, respectively. 
This feature remains true for all physical VEVs and shows that the only physical particle 
in the Schwinger model is the massive Schwinger boson. 

Further it may be seen easily that all VEVs fulfill the vacuum clustering property, e.g. 

lim (S ± (x)S ± (y)) = (S ± (x)) (S ± (y)) . (43) 

x—y—>oo 

In fact, vacuum clustering is another reason that makes it necessary to introduce an instanton 
vacuum. 

The above calculations may be easily generalized to higher VEVs of chiral densities. The 
general formula reads (see e.g. fl3"Tf ) 



(S Hl {xi) ■ ■■S Hn {x n ))o = e ifce (-) ri exp[^a i a j 4vrL) / , (xi - Xj)] (44) 

i<j 



k = Oi = n + — n_ 
i=i 

where cij = ±1 for Hi = ±. This result we need for the mass perturbation theory. 
Further VEVs that may be easily found are the field strength VEV 
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(F(x)F(y)) = ^D x D y (/3(x)/3(y)) = -^D m {x - y) - 5(x - y) (45) 

and the vector current correlator. For the latter it is most efficient to introduce a vector 
current source into the path integral, 

e£ M = e^d v b (46) 



eB.J 11 = -^e^be^dx® = =: <f>A (47) 

V 71 " V 71 " 

where A = -4=D6 is the source for the Schwinger boson $. Now the inclusion of the source 

V 71 " 

consists in the substitution (3 — > (3 + b in the second (anomaly) term of the effective action 
(36) and in the exact fermion propagators (28) and zero modes (30). The remaining path 
integral computation is similar to the one which we discussed previously and leads to the 
following VEVs (for n chiral densities, which are at the same time generating functionals 
for the Schwinger boson) 

(S Hl (xi) ■ ■ • 3 Hn (x n )) [X] = e ifc9 (-) n exp[^a l a i 47rD /t0 (x i - Xj ) ■ 

i<j 



exp 



- / dy 1 dy 2 X(y 1 )D IM) (y 1 - y 2 )A(y 2 ) + 2i^J2^i I dy\{y)D tXo {y - (48) 
J i=i J 



(for an explicit computation see p3| ). 

Observe that, although b was needed in the intermediate computations, the final result 
can be expressed entirely in terms of A, and only massive propagators occur. This again 
shows that the massive Schwinger boson is the only physical particle of the theory. 

After this short discussion and formulae collection of the massless Schwinger model we 
are prepared for the mass perturbation expansion of the massive model. This will be done 
in the forthcoming sections. 



V. VACUUM FUNCTIONAL AND VACUUM ENERGY DENSITY 

By simply expanding the mass term in (20), e™/'™, the Euclidean vacuum functional 
for the massive Schwinger model may be traced back to VEVs of the massless model and 
some space-time integrations ( fl27|), 

oo n „ n 

Z(m } 9) = Y, —r / dx x . . . dx n {\{ ^(x,)^(x 4 ))o (49) 

where (117=1 S(xi))o = ([\2=i(S+(xi) + S , _(xj)))o may be inferred from (44). There, all 
contractions among SuXx^SH^Xj) produce exponentials of the massive propagator _D M0 (xj — 
Xj). For a first, rough approximation we may use the fact that the massive scalar propagator 
D^x) vanishes exponentially for large argument. Therefore, when integrating over space 
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time and expanding the exponential, all contributions from D 1 will be ignored for the 
moment, supposing that the space time volume V is sufficiently large. In this case the 
integrations in (49) just produce factors of V. Further, when inserting (44) into (49) we 
have to sum over all possible distributions of n + = n — n_ pluses and n_ minuses on n scalar 
densities S. This results in a factor ) . Therefore we find for the n-th order term 



i=i 



Yl / dxiS(xi)} 



n_=0 



n 



j(n— 2n_ 



(EFcosi 



(50) 



and for the normalized vacuum functional 



Z(m, 



2(0,0) 



exp(m£Vcos( 



(51) 



(we ignored terms like m n V n ~ l in this approximation compared to m n ~ 1 V n ~ 1 , therefore (51) 
is the first order result in m). This result is wellknown, and its consequences for the vacuum 
structure and spectrum of the Dirac operator are discussed in great detail in [[4^ . 

As we have seen, the exponentials e ±47rZ)Mo ^ produce volume factors V upon integration, 
and, consequently, they will not lead to an IR-finite perturbation expansion for V — > oo. 
Therefore, we have to expand the (products of) exponentials in (49) into the functions E±(x) 

E ± (x) := e ±4 " D "oM _ l 



E^\x) := e ±4wD "o(-) 



n i 

Et[(±4^ () (x)) Z 
i=o L 



E ( ±\p) = J d 2 xe ipx E { ±\x) , E± := / dxE ± (x) (52) 

(where we defined some related functions for later convenience). The E±(x) decay exponen- 
tially for large \x\. 

Inserting the exact VEVs (44) into (49) and using the notation (52), we obtain for the 
vacuum functional Z(m, 9), order by order 



n=l: 



n=2: 



^Jdx(e^ + e^) = ^V(e^ + e^) 



(53) 



m 2 / E \ 2 



2! V2 



(-) / dx x dx 2 ^O^D^-xi) + 2e -^D^{ X1 -x 2 ) + e -2ie e 4nD^ Xl -X2) 



m 2 / E\2 



2! V2 



(±) z [^2 (e 2^ + 2 + e -2 i6) + V ( E+e 2ie + 2E _ + E+e -xe) 



(54) 



n=3: 
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m 3 / E \ 3 



V 2 (3E + e 3ie + 3(£+ + 2E_)e ie + 3(£+ + 2E_)e~ ie + 3£ +e - 3ie ) + 

\/((3£ 2 + E + x £ + x £ + )(e 3if? + e~ 3ie ) + 3(2£+£_ + ^+ £ + xlx E_)(e ie + e - ^)) 

(55) 

n=4: 

m 4 / S \ 4 



4! V2 



-2«6» | „-4i6>\ 



\/ 3 (6E + e 4ie + 12(£+ + £_)e 2lf? + 12(£+ + 2£_) + 12(E+ + EJ)e~ 2td + QE + e~ m ) + 
V 2 ((15E 2 + + 4£ + x E + x £ + )(e 4if? + e~ 4ie ) + 



4(3£^ + 9E + E_ + 3E! + E + x E + x E + + 3E + x E_ x E_)(e 2i6 + e~ 2id ) + 



+ Q(El + 8E + E_ + E 2 + AE + x x + ... 



(56) 



where the cross indicates convolutions, e.g. 

E+x E + x E + = J dy 1 dy 2 E + (y 1 )E + (y 1 + y 2 )E + (y 2 ) 



(57) 



and we displayed the result up to the accuracy we need. Observe that the result is not 
obtained by just expanding polynomials like (1 + E + (xi)) 1 , because e.g. a third power in 
E + (xi) may contribute to V n ~ 3 E\ or to V n ~ 2 E + x E + x E + . Concerning the dimensions, 
observe that E±(x) is dimensionless and, therefore, e.g. [E±] ~ [V], [E± x E± x E±] ~ [V 2 ], 
etc. 

In a next step we rearrange the terms (53) - (56) in rising powers of V: 



^h|(e^ + e~*) + ^y(E + e™ + 2E_ + E + e~^)+ 



m 3 / E \ 3 



— (-) ((3E 2 + +E + xE + xE + ) (e m +e- 3%e )+3(2E + E^ +E 2 _ +E + xE_xE_) (e ie +e' w )) +. . . 



+ 



V 2 r 2 /S x 2 
-2\ m ' 



^\ e ™ + 2 + e~ 2W ) + m 3 &(E + (e m + e~ m ) + (E + + 2E„){e w + e' w )) + 



V2 
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m 4 f E \ 4 



12 ^2 



(-) ((15£ 2 + 4E + xE + x E + )(e AW + e~ m ) + 



4(3£ 2 + 9£+£_ + 3£ 2 + E + x E + x E + + 3E + x x £_)(e 2ie + e" 2if? ) + 



6(£ 2 + 8£ + £_ + E 2 _ + 4£+ x £_ x £_)) + . . . 



+ 



4 -Ex4 



rn 

~2 V2 



(-) (3£ + (e 4 * e + e~ 4 ^) + 6(£+ + E_){e™ + e" 2 *) + (E + + 2£_)) + 



+ 



4! 



+ . . . 



V V 2 2 V 3 3 

This result indicates an exponentiation of the exact vacuum functional, too, 

Z(m,9) 



Z(0,0) 



3 Ve(m,«) 



where 



S T7i 2 S 2 

e(m,6) = m-2cos0 + — (-) (2£+ cos 2^ + 2£_) + 



(58) 



(59) 



— (-) ((3£ 2 + x x £ + )2cos3# + 3(2£+£_ + E 2 _ + E+ x E- x £_)2cos#) + . . . 

(60) 

This exponentiation, equ. (59), is, of course, very important, because it guarantees that 
physical VEVs do not depend on the space-time volume V and are therefore IR-fmite. 
Therefore, it would be nice if equ. (59) could be proven generally. This is indeed possible 
and shall be discussed next. 

For this purpose, let us write Z(m, 6) as a double sum, 



oo I I 

zfo*) = E ir E ^,». 

1=0 L n=0 

where the factorial has been introduced as in the perturbation expansion (49). 
Because of the multinomial formula 



(61) 
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m 



(e(m,e)r = (mc hl + —c 2A + ..T 



the exponentiation condition reads (the ^7 stems from the factor ^) 

(n + Z)! ~ d "V(n + Z)! 



or 



Cn+Z,n 



■A (n + /)! /Ci ; i\fel / \K+l , , 





Next we have to answer the question where the volume factors V come from. The first power 
V stems from the fact that all "propagators" E±(xt —Xj) depend on coordinate differences. 
Higher powers occur when the corresponding "Feynman graph" is disconnected. Let us show 
a graphical example in Fig. 1. 



vAAA/ 

Fig. 1 



Here each wavy line represents a E±(x) and their endpoints represent "vertices" 2 rpe ±i6 '. 
These Feynman rules will be discussed in detail in Section 7. 

Obviously, Fig. 1 is of third power in V and of 8th power in m (8 vertices). In general, 
ci. n is just the sum of all graphs of Z-th order that consist of n connected pieces. 

Therefore, we just have to prove that (63) is the correct combinatorial formula for the 
distribution of n connected graphs of total order n + I on n + I vertices. 

But this is easy to understand. Consider e.g. a graph q 1+ / 2i 2 consisting of two connected 
pieces Q l5 i, q 2i i. There are (h + l 2 )\ possibilities to distribute the two connected graphs on 
h + h vertices. However, l ± \ (l 2 \) ways exist to rearrange (q 2j i) on a given subset of l\ 
(l 2 ) vertices, therefore these factors must be divided out. This leads to 

Ch+i 2 ,2 = (h + hV- y 

and easily may be generalized to explain (63) up to the ki factors. When some connected 
subgraphs c itl occur more than once, ki > 1, there is an additional symmetry factor that 
counts for the fact that an exchange of identical subgraphs c it i leads to the same contribution 
to c n+ i tTl . This explains formula (63). 
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So we have proven the exponentiation of the vacuum functional, (59) and, as a conse- 
quence, the IR-fmiteness of the mass perturbation theory. In fact, the expansion of e ±47rZ V> ( x ) -. 
(1 + E±(x)) into E±(x) is analogous to the cluster expansion of statistical physics models. 
There, our result correspond to the fact that the free energy is an extensive quantity. A 
more rigorous discussion of these features can be found in p^j30|1 . 

Before ending this section we want to give an explicit (numerical) formula for the lowest 
orders of the vacuum energy density, e(m,9), and, as a consequence, we will meet the 
problem of UV regularization. For a numerical evaluation of order m 2 we have to compute 
the coefficients E< and E_. First, both E + and E_ are proportional to \. Scaling u out, 

Mo 

we find 

f4E + = fi*[ d 2 xE + (x) = J d 2 x(e- 2Ko ^ - 1) 

/■oo 

= 2tt / drr(e- 2Ko{r) - 1) = -8.9139 (64) 
Jo 

E + (x) is well behaving (E + (0) = —1), so the numerical integration is straight forward. 
E-{x) is singular at x — 0, E-(x) ~ ^ for i -> 0, but this singularity can easily be 
understood and removed in a unique way. Indeed, this singularity is just the free fermion 
singularity, as can be seen by rewriting E-(x) like (see e.g. equ. (41)) 



(f) V (x) + 1) = (S + {z)S.(0)) = « Gl(x) = ^- 2 . (65) 



This singularity may be isolated by a partial integration: 

dr 



vlE_ = [ d 2 x(e 2Ko ^ - 1) = 2tt f°° ^(e 2 ^M+2inr _ ^ 
J Jo r 

= 27r[(lnr)(e 2 ^W +21nr -r 2 )]" o + 

2tt / dr2(\nr)((K l (r) - _) e 2 ^M+2inr + s ^ 
Jo r 

(K' q = —Kx). Observe that the first term precisely leads to the free fermionic singularity at 
the lower boundary (and vanishes at the upper boundary). So the second term is the unique 
and finite result we are looking for. The numerical integration gives 

fi 2 E_ = 9.7384 (67) 

In the literature there exist other UV regularizations as, e.g. the introduction of an 
additional Thirring-type interaction ( []30|j3~3]0 . Of course, when this Thirring-type coupling 
constant is set equal to zero at the end of the computations, the result agrees with our 
regularization. 

With the help of (64), (67) the vacuum energy density reads 

e(m,^) = mScos^ + -^-(-) (-8.9139 cos 26 + 9.7384) + o(m 3 ). (68) 

Higher order contributions to e(m, 8) (and to bosonic n-point functions to be discussed in 
later sections) may contain logarithmic singularities like in (66) and may be regulated in a 
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similar manner. However, one may perhaps hope for an even better UV behaviour in higher 
order computations. After all, the model is known to be super-renormalizable (in ordinary 
electrical coupling perturbation theory this feature is obvious, as the vacuum polarization 
graph (= a loop of two fermions) is the only divergent diagram). So one might expect 
cancellations of divergencies in higher order computations. And, indeed, for the o(m 2 ) 
contribution to the Schwinger mass precisely this happens (see Section 8, (99)). A more 
detailed discussion of this problem will be given elsewhere. 



VI. CONDENSATES, SUSCEPTIBILITIES AND CONFINEMENT 

In this section we will try to extract some physical information from our results obtained 
so far. The simplest task is the computation of the fermion condensate, 

< SW > = <*(,)*(,» = L^A-zw) = A e ( )M) 

— TTl 

(V(x)ty(x)) = Yl cos 9 + — Y?(E + cos 29 + E_) + o(m 2 ). (69) 

Of course, the order zero result is the condensate of the massless Schwinger model. From 
this result the pseudoscalar VEV (P(x)) may be obtained, e.g. by the substitution cos n9 — > 
i sin 710 (n — 0, 1, . . .), 

— in 
(P(x)) = (*(ar)7 5 W(a;)) = i£ sin 9 + i— £ 2 £+ sin 29 + o(m 2 ). (70) 

Another quantity that may be easily obtained is the field strength condensate (see equ. (21)) 

2tT U 2"7T 7T 

I Fix)) = — ^(m, 9) = roS sin 9 - -m 2 Y?E+ sin 29 + o(m 3 ). (71) 

e o9 e e 

Therefore, as discussed in Section 2, whenever there is a classical background field, F ~ 9, 
there remains some effect on the quantum level and the screening is not complete in the 
massive model. 

Observe that there is a relation between (P) and (F). This is due to the anomaly 
equation (11), 

= d li {Jg) = -{F)-2im{P) (72) 

7T 

(one-point functions are always constants because of translational invariance). 

By performing a second derivative on e(m, 9) one may infer the susceptibilities. For the 
scalar susceptibility we get 

Xs = J dx{S(x)S(0)) c = ^e(m,9) = \^{E + cos 29 + E_) + o(m) (73) 

(the subscript c indicates the connected component; our discussion of the last section ex- 
plains why only connected components may contribute to e(m, 9)). For the pseudo-scalar 
susceptibility the mixed contractions S+, S- get a minus sign, 
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X P = J dx(P{x)P{0)) c = ~£ 2 (£+ cos 26 - E.) + o{m). (74) 
The topological susceptibility is 

Xtop = (=^f / dx(F(x)F(0)) e = ~<m,B) = 



mE cos 6 + m 2 Y?E + cos 20 + o(m 3 ). (75) 

There is a relation between x p and Xto P ; however, this and related questions will be discussed 
in a later section, where the Dyson- Schwinger equations corresponding to the equations of 
motion (10), (11) are derived. More on the physical meaning of susceptibilities in two- 
dimensional models may be found e.g. in [|43| . 



A further physical feature that we are able to discuss by using our results obtained so far 
is the confinement behaviour of the massive Schwinger model. This was recently discussed 
in in first order mass perturbation, and we follow their approach and generalize their 
result to arbitrary order. 

A usual way to derive the confinement property is the computation of the string tension 
from the Wilson loop. The Wilson loop for a test particle of arbitrary charge g = qe is defined 
as (the additional factor i in Stokes' law is due to our Euclidean conventions (imaginary F), 



see e.g. [IE, 1 1 



W D = (e i9 ^9D A ^) = ( e 9j D F ( x ) d2x ) = ( e ^iqf D v(x)d 2 x^ ^ 

where v(x) is the Pontryagin index density. Further dD is the contour of a closed loop and 
D the enclosed region of space-time. We are interested in the string tension for very large 
distances; further we are able to explicitly separate the area dependence, therefore we may 
set D — > V in the sequel. 

For the VEV of an exponential the following formula holds, 

(e— ) = ex P [E^ ! K) c ] (77) 

n=l n - 

where () c denotes the connected part of the n-point function. These VEVs are given by 

r d n 
(u n ) c = Vjdx 2 ... dx n (v(0)u(x 2 ) . . . u(x n )) c = V(-i) n —e(m, 6) (78) 

as is obvious from the vacuum functional (21) (as usual, disconnected VEVs are of higher 
order in V). The vacuum energy density of arbitrary order may be written like 

oo 

e(m,6) = ^e*cos/# (79) 

1=0 

where we used the symmetry of e(m, 9) with respect to positive and negative instanton num- 
ber, and ei acquires contributions from instanton sectors k = ±1. Performing the derivatives 
(78) we have to separate even (~ cos/0) and odd (~ sinW) powers of derivatives. We find 
for the Wilson loop 
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w 



2n— 1 oo 



n=l 
oo 



(2n)! ^ } „tl(2n-l)! 



^(_l)"/2n-l ejSinW 



2=0 

oo 



2=0 



exp[vE^os^E(-ir^P + 

2=0 n=l \ Z "'J- 2=2 

oo oo 

exp [V e/ cos W(cos 2ixql — 1) — V ^ e/ sin /# sin 27rgZ 
2=0 /=o 



00 00 (2iral) 2n ~ 1 - 

S (2n-l)!- 



(80) 



The string tension is defined as 



(7 



J oo 

— In = ^2 e 2 ( cos ^(1 — cos 27rg/) + sin 19 sin 27rg/ 

^ 2=0 



(81) 



and may be interpreted as the force between two widely separated probe charges g = qe, 
where all quantum effects are included. 

We find that, whenever the probe charge is an integer multiple of the fundamental charge, 
q G N, the charges are screened, and the Wilson loop does not obey an area law. Observe 
that this result is exact ! 

For noninteger probe charges there is no complete screening and the string tension may 
be computed perturbatively, 



a = mTj(cos6(l — cos27rg) + sin#sin27rg) + o(m 2 



(82) 



showing that in the massive Schwinger model and for noninteger probe charges there re- 
mains a constant force (linearly rising potential) for very large distances. So in the massive 
Schwinger model true confinement is realized instead of charge screening in the general case. 
Observe that this formation of a long range force is a strictly nonperturbative phenomenon 
in the sense of conventional perturbation theory, because only nontrivial instanton sectors 
(k ^ 0) contribute to the string tension, as is obvious from (78). 

On the other hand, in the massless model arbitrary probe charges are screened. (Discus- 
sions of screening and confinement in the massive Schwinger model within other approaches 
may be found e.g. in |TTi^4|,^8i^0|,^5i , and the behaviour found there agrees with our result.) 



VII. FEYNMAN RULES FOR MASS PERTURBATION THEORY 

The interaction Lagrangian for the mass perturbation expansion reads L\ = m\I/\I/, see 
(49). On the other hand, the formulae (44), (48) for VEVs of the massless model, that we 
need for the perturbation expansion, contain chiral currents S± instead of the scalar one in 
L\. As a consequence, each vertex corresponding to L\ contains in fact two vertices, 

m(S(x)) = m(S+(x)) + m(S-(x)) = ^e w + ^e~ ie (83) 

corresponding to the two chiralities. 

Further, these two types of vertices are connected by two types of propagators, namely 
S + (x)S + (y) and S-(x)S-(y) by E + (x — y), and S + (x)S_(y) and S-(x)S + (y) by E_[x — 
y). Further, because all vertices may be connected to each other (see (44)), up to n — 1 
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lines E±(x — yi) may run from one vertex x to the other vertices yi for a n-th order mass 
perturbation contribution. 

As a consequence, the Feynman rules acquire a matrix structure. More precisely, the 
propagator corresponding to the E±{x) is a matrix, which in momentum space reads 

e (p) = (E + (p) E.(p)\ 

{P) \E_{p)E + {p)) {M) 

where the individual entries correspond to the individual (Si(x)Sj(y))o, i,j = ±, propaga- 
tors. 

Each vertex, where n propagator lines £{pi) meet, is a n-th rank tensor Q . Only two 
components of this tensor are nonzero, namely 

&0++...+ = — e , C, = — e (85) 

(corresponding to S — S + + S-). E.g. the vertex where two propagators meet is a matrix 

(mS id n \ 
2 o (86) 

Internal lines must be £{p) propagators; external lines, however, may be boson lines, too, 
when we treat bosonic n-point functions (i$(xi) . . . i$(x n )) m . In (48) we see that each boson 
that is connected to a S- vertex acquires a minus sign. Therefore, the rule for a boson line 
is that each vertex S = S + + S_, where a boson line meets, is multiplied by times the 
boson propagator p2 ~^ times the pseudoscalar vector P, where 

When n bosons meet at one vertex, one may, instead of contracting that vertex with n 
vectors P, contract it with one P (S) if n is odd (even). Of course, the number of indices 
of the vertex must be reduced accordingly. 

These Feynman rules may be given by the graphs of Fig. 2 (we display them in momen- 
tum space). 

D M ... 

Go 



vAAAAA/ 



i x 

Fig. 2 

Here Q denotes the renormalized coupling that may be found like follows. The bare 
couplings are m{S±)o = H ^e ± * e , therefore the renormalized couplings are defined as 
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go = m(S + (x)) m , g* e = m{S-(x)) m 



(88) 



and Q is constructed out of go, g% like Q out of H ^e ±l61 , see (85). Graphically, Q is just the 
sum of all graphs that may be attached to the bare vertex Q Q , see Fig. 3, 



. = . + t + r + 



Fig. 3 

i.e. the sum of all graphs where a line of propagators either starts and ends at the same 
point, or it ends in the vacuum with zero momentum, go may be easily computed to be 

ge = ^ + (^)\e™E + + E_) + o{rn*) 

=■ 9i + 92+ o(m 3 ) (89) 
In fact, because of go = m{S + ) m , it is related to the vacuum energy density, see (69), 

d 

9e + 9*e = m—e(m,9). (90) 

Further, because each vertex may be connected to all the other vertices in our theory, the 
contributions to the renormalized coupling, Fig. 3, may be attached to each vertex of an 
arbitrary graph and, therefore, all the vertices of the theory may be renormalized from Q 

to g. 

As a further example of these Feynman rules we will investigate the bosonic two-point 
function in the next section (we will compute it explicitly there). In addition, we will use 
these Feynman rules extensively in Sections 9-11. 



VIII. THE SCHWINGER MASS 

As usual, in order to compute VEVs for the massive Schwinger model one has to insert 
the corresponding operators into the path integral (20) and divide by the vacuum functional 
Z(m,9) (see (59), (60)), 

1 00 m n n r 

(0)m = yT—rfAO £ — [J / dxMxiMxi)^ (91) 
Z(m,U) n=0 n\ i=1 J 

We will find that via the normalization all volume factors cancel completely, as it certainly 
has to be. 

The Schwinger mass may be inferred from the Schwinger-boson two-point function. There- 
fore, as a starting point for the mass perturbation computation we need the following n-point 
functions of the massless model, 

n y 

(i*(ife)i*(yi) n^te^o = e lfce (-)Y^<^ D ^-0. 

,■—1 A 



23 



D, ( yi -y 2 ) + 4n(J2H ai D«>(xi-y2))(T,(-) ajD »o(xj-yi 

3=1 



i=l 



(92) 



i=i 



which may be easily computed from the generating functionals (48). 

For a perturbative computation of the Schwinger boson propagator we simply have to 
insert successive orders of (92) into the perturbation formula (91). Doing so, we find up to 
second order 



(i$(2/i)i$(j/ 2 ))r 



1 



Z(m,0) 



D, ( Vl - y 2 ) + m-(e w + e'^VD^y, - y 2 ) + 



A7rm-(e id + e' ie ) J dxD^x - Vl )D^{x - y 2 ) + 



-^-(-) {e 2ie + e~ 2i9 ) J dxidx^D^ - y 2 ) + 4tt( j D w (x 1 - Vl ) + D^{x 2 - Vl ))- 



■{D^x, - y 2 ) + D^(x 2 - y 2 ))]e 4 ^o(« 2 ) + 



^j-(^) 2 2 J dx 1 dx 2 [D /t0 (y 1 - y 2 ) + An(D^(xi ~ Vi) ~ D m {x 2 - y x ))- 



{D^{ Xl - y 2 ) - D^{x 2 - y2 ))] e -^ ( xl - X2 ) 



(93) 



Inserting Z(m, 6) up to second order (see (59), (60)) and expanding the denominator in the 
usual perturbative fashion, we arrive at 

(i^(y 1 )i^(y 2 )) m = D^ ( yi - y 2 ) + 47rm£cos# J dxD^ (x - yi)D w (x - y 2 ) + 



/ S 2 f 

4?rm 2 (-) cos 26 J dx 1 dx 2 [2D )l() (x 1 - yi )D^( x i - V2) + 2-D Mo (xi - y 1 )D /Xo (x 2 - y 2 )\- 



X) 2 [' 

(E + ( Xl - x 2 ) + 1) - 87rm 2 (-) cos 20 V J dxD m (x - y^D^x - y 2 ) + 



Yl 2 f 

4vrm 2 (-) J dxidx^D^fa - yi)D fl0 (x l - y 2 ) - 2D fI0 (x 1 - y^D^fa - y 2 )\- 



Yl 2 t 

{E-{x x - x 2 ) + 1) - 87rm 2 (— ) V I dxD^ix - yi )D^(x - y 2 ) 
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2 f 

= D^iyi -y 2 ) + 47rm- cos6> I dxD^(x - y^D^x - y 2 )+ 

£ 2 f 

87rm 2 (— ) J dx 1 dx 2 D fl0 (x l - j/i)D /J0 (x 2 - y 2 )[cos26' J B + (xi - x 2 ) - - x 2 )]+ 

S 2 /" 

87rm 2 (-) (cos20£+ + £_) y dxD^x - y^D^x - y 2 )+ 

S 2 /" 

87rm 2 (— ) (cos26> - 1) / dx 1 dx 2 D fl0 (x 1 - y^D^fa - y 2 ) (94) 

where E±,E±(x) are given in (52) and we used the x — > —a; symmetry of all occurring 
functions. The last term stems from a disconnected part of (93) and must be subtracted. 
Observe that, as claimed, all volume factors V have dropped out. 

We may easily check that (94) is the right expression by depicting the corresponding 
Feynman diagrams in Fig. 4 




+ 471 ( ^ ^ + ' x + 

— I — + + ... ) 

Fig. 4 

and find that (94) and Fig. 4 coincide. Again, the last graph is disconnected and must 
be subtracted. 

In order to obtain the second order result for the Schwinger mass, we rewrite expression 
(94) in momentum space and substitute all functions by their Fourier transforms (thereby 
the convolutions turn into products; the disconnected term in (94) is omitted), 

— c -1 1 
(i$i$) m (p) = — 2 +47rmScos6»- 



P 2 + ti (p 2 + ^) 2 



27rm 2 £ 2 1 [cos2fl(E + + E + (p)) + E. - E. 

(P + fJ, Q ) 



j-^fl -47rmScosg — 9 - 27rm 2 £ 2 [cos 26(E + + E+(p)) + E- - EJp)] . ) 
P^ + K K P 2 + K P 2 + K J 



P 2 + Vl + 47rm£cosfl + 27rm 2 £ 2 [cos2#(£ + + E+(p)) + E_ - E(p)\ + (47rm£cos#) 2 ^^ 
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+ o(nf) (95) 

Therefore, for finding the mass pole, p 2 has to obey the equation (after a rescaling p' 2 = 
E' ± = E ± ( l x 2 = l) = f x 2 E ± etc.) 

,2 TnYl „ m 2 £ 2 „ r . ~. . .. 47T . 

p' = -1 -4tt— cos#-2tt — _cos20[E+ + ^ + (p') + -^ r]- 

/i p 11 p + 1 

m 2 E 2 ~ , „ 4tt 



2tt 



— [E'_ - E_(p') + - T —]. (96) 

The second order part (the term in square brackets) may be further evaluated like (for the 
cos 29 term) 

[•••] = / d 2 x[e~ 2KM - 1 + ^(e-atfoCM) _ i + 2K (\x\)} = 

roo p2,tv 

/ drr[27r{e- 2Koir) - 1) + / dfle i|p l"°^( e -2^o(r) _ x + 2 K (r))} = 
Jo Jo 

2tt / drr[e- 2Ko{r) - 1 + J {\p'\r){e- 2Ko{r) - 1 + 2# (r))] (97) 
Jo 

where Jo is the Bessel function of the first kind. This expression behaves well around \p'\ = i 
and therefore we may set \p'\ = i because deviations from this value are of higher order in 
m. Using Io( r ) — Jo{^ r ) we find 

— [...]=: A = / drr[e- 2Kai - r) -l + h{r){e~ 2Ka{ - r) -l + 2K {r))} = -0.6599 (98) 
2n Jo ' 

Analogously we find for the other second order term (containing the EJ) 

_[. ..]=:B= / drr[e +2Ka{r) - 1 + I {r){-e +2Ko ^ + 1 + 2K (r))] = 1.7277 (99) 
2ix Jo 

In this expression (99) the nice feature of cancellation of UV divergencies occurs. Indeed, 
both e 2K °^ and Io(r)e 2K °( r > diverge like \ for small r (this divergency corresponds to the 
free fermion field divergency of the underlying theory that we discussed in Section 5, see 
(65)), but obviously the divergencies cancel each other. In fact, this cancellation was already 



observed twenty years ago in [[T(J within a bosonization approach. 



Collecting all results we find for the Schwinger mass in second order 



/ = 4 = 1 + 4tt— - costf + A^(-)\a C os29 + B) (100) 
Ho hq hq Ho V /V 



or, inserting all numbers (remember — = equ. (40)) 



m . m 2 m 2 



^ = ^(l + 3.5621 • — cos 9 + 5.4807 ■ — - 2.0933 ■ — cos 29). (101) 
A'o f*o Ho 



For the special case 9 = our result (101) precisely coincides with the result in | 4"6"|| , where 
the second order correction for 9 = was computed within bosonization and using near light 
cone coordinates. In the same article this result was compared to a lattice calculation ( [|52|), 
and a good agreement is obtained within the range of the expansion parameter ^ where 
the lattice calculations were performed. 
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IX. DYSON-SCHWINGER EQUATIONS AND EXACT iV-POINT FUNCTIONS 



The Schwinger boson - which we discussed in the last section - is an interacting particle 
in the massive Schwinger model. As a consequence, we will find that it forms n-boson 
bound states. In principle, their masses could be computed analogously to the previous 
section, by a computation of the corresponding Schwinger-boson 2n-point functions and by 
the determination of their mass poles. 

Here we will adapt a slightly different method. We will discuss the Dyson-Schwinger 
equations that follow from the equations of motion (10), (11), acquiring thereby a deeper 
insight into the structure of the model. With the help of these Dyson-Schwinger equations 
we will be able to re-express the n-point functions of the model in a way that is more suitable 
for our discussion. We will find in this way that the spectrum of the theory is even richer 
than expected. There is a second stable particle in the theory in addition to the Schwinger 
boson, namely the two-boson bound state, and unstable higher bound states may be formed 
out of both these stable particles. Further, we will find that both particles may occur in 



final states of decays and scattering processes (see also |62|-|64[). 

When we use the generating functional for Schwinger bosons (48) for a computation of 
the Schwinger boson 2n-point function in lowest order in m, we find a contribution that may 
be depicted graphically like in Fig. 5, 



Fig. 5 

and the vertex corresponds to a coupling constant c, where 

c = (4vr) n mS cos 9. (102) 

Therefore, this lowest order coupling mediates an attractive force for \9\ < |, and in this 
range of the vacuum angle 9 the formation of bound states has to be expected, at least for 
sufficiently small fermion mass m. The restriction \9\ < £ shall be assumed in the sequel. 

In the introduction we wrote down the two equations of motion that relate the field 
strength operator and the fermionic vector current operator, namely the Maxwell equa- 
tion (10) and the anomaly equation (11). By introducing the Schwinger boson (9) and by 
eliminating the field strength we may derive the equation (expressed in real fields P, i$) 

M x i$(x) := (n x - /4)i<£>(x) = 2^mP(x). (103) 

(where we introduced the operator M x for convenience). 

[Remark there is a slight difference between the equations (10), (103) and equation (11). 
Whereas the anomaly equation (11) holds as an operator relation, and consequently on all 
states, (10) and (103) are only true on physical states (i.e. on states that are invariant with 
respect to small and large gauge transformations; for a deeper discussion of this feature see 
e.g. P,[TT|1). However, the introduction of the 9 vacuum within our path integral approach 



precisely corresponds to the introduction of the physical vacuum (see Section 2). Therefore, 
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for our physical VEVs we may use equ. (103), as we will verify by explicit perturbative 
computations.] 

Next we define the amputated, connected bosonic n-point functions (FT. . . Fourier trans- 
form) 

M W (pi, ..., Pn ) = FT(M X1 . . . M Xn {i<S>( Xl ) . ..i$(x n )) c J. (104) 



A. Two-point function 



So let us study e.g. the connected two-point function (i&(x 1 )i$(x 2 ))m- Graphically it 
may be depicted like in Fig. 6. 




+ — + — * + ... ) 

Fig. 6 

We find the following behaviour: all graphs where both boson lines meet at one and the 
same vertex contain just the corrections that change this vertex from the bare one Go to the 
exact (renormalized) vertex Q, see Fig. 3. In all the other graphs, each individual vertex is 
renormalized in the same way, too. 

So let us re-express Fig. 6 in terms of the renormalized coupling, and, in addition, 
amputate the two external boson lines (and the pseudoscalar vectors P, see (87)). We 
obtain Fig. 7, where we introduced the exact propagator that is defined in Fig. 8. 



-f \/^) /~ = ~f /- + 471 *wwvQaamm 

Fig. 7 

»wwvQaaaw« — • + •wwv* +*ww*ww^« -|- <^wv\i» -|- 



-)- l^^^p -|- rfj5w3>Sji -|- rfvw^J3•^AAAA<• -|- 



Fig. 8 

Here it is understood that the left and right vertices of each graph in Fig. 8 are the initial 
and final ones where we amputated the bosons. We introduce for the above exact propagator 
of Fig. 8 the name GH(p) in momentum space (matrix multiplication is understood, see the 
Feynman rules of Fig. 2), because we will need it frequently 

gu( P ) -.= g + gs( P )g + gs( P )gs( P )g + ... (105) 



28 



This exact propagator H(p) consists of a constant part (the sole vertex in Fig. 8 without 
£(p) lines) that is a scalar (because two boson lines meet on this one vertex), and of a true 
two-point function (depending on p) , where the initial and final vertices are pseudoscalars. 
Inserting all factors properly, Fig. 7 may be written like (M^{p,p) = M< 2 >(p)) 

M {2 \p) = -(p 2 + 1$) + 4vrm( 1 S) m + Anm 2 (Ppfjp) (106) 

where 

m(S) m = P T gP = S t g t =g e + g* e (107) 

m 2 (PP)Jp) ee P T g(U(p) - 1)P (108) 

where P T = (1,-1) is the transpose of the vector P, (87), and matrix multiplication is 
understood (the single vertex Q we may interpret either as a two-component object that is 
contracted by two vectors P or as a one-component object that is contracted by one vector 
S). 

This is just the momentum space version of the Dyson-Schwinger equation for the two- 
point function, 

M yi M y2 (i<Z>( yi )i<f>(y 2 )) = M yi 5( yi -y 2 ) + 

A7im(S( yi )) m 5( yi - y 2 ) + Aixm 2 {P{ yi )P{ y2 )Y m . (109) 

Observe that the (P(yi)P(y2)) propagator includes an arbitrary number of bosons propa- 
gating from 2/i to y 2 , even in least order. This will be important in the sequel. 

The key observation for the computation of bound states is the fact that the exact 
propagator H(p) may be resummed. This resummation relies on the following observation. 
All diagrams that fall into two pieces when they are cut at a vertex, factorize in momentum 
space, i.e. they are a product of two functions of p, see Fig. 8. The opposite type graphs 
are called non-factorizable (n.f.). 

Here we should be more precise about the cutting. We stated in Section 7 that the 
vertices are tensors, so how to cut such a vertex? Suppose e.g. we have a vertex where three 
lines meet and we want to cut it in a way that two wavy lines belong to the left hand side, 
and one line to the right hand side. Then we rewrite the vertex like 

Gijk = SijiGii'Si'k , i,j,k,l,V = ± (110) 

where Qw is the vertex matrix (85,88) and the 5^...^ are generalizations of the Kronecker 
delta Sij, i.e. 

S++-+ = 5 = 1 , 5 i:L ...i n = otherwise (HI) 

Therefore, we may write for the sum of non-factorizable graphs, which we call A (see Fig. 
9) 
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Fig. 9 

/cPq 
j2^8 ik k'£ki(q)Gw£i'm(q)£k'm'(q-p)8jmm' + ■■■ (112) 

The matrix A(p) may be rewritten like 

A(V) = ( ^±?+)n.f.(P) ( S ±^-)n.i.(P) \ ( U3 ) 

\(ss + ) q Ap) (ss-) q Ap)J 

The entries of this matrix are, however, related (e.g. (SS-) nf (g$,p) = (S + S + ) n f {.9e,p)-> 
as may be checked from the perturbative expansion) and, therefore, we find for the product 
QA{p) (which we need in the sequel) 



®(ge,p) = ge(S+S + ) nf (g e ,p) , /3(g g ,p) = ge(S + S-) nf (g e ,p). (115) 

Now we may collect all n.f. graphs in (105), Fig. 8, e.g. on the left hand side, and find that 
they are again multiplied by all graphs that occur in Fig. 8. Therefore we may write for 
gU(p) of equ. (105) 

gu(p) = g(i + A(p)gu(p)). (116) 

Equation (116) may be solved for Tl(p) by a matrix inversion and has the solution 

TIM - — ( 1 " a ^ eiP) ^ 9 * 9,P) ) (117) 

where N(p) is the determinant of the matrix that had to be inverted, 

N(p) = det(l - gA(p)) = 1 - a(g e ,p) - a(g* e ,p) + a(g e ,p)a(g* e ,p) - (3(ge,p)f3(g* e ,p). (118) 

We will find that the zeros of the real part of the denominator N(p) will give us all bound- 
state masses, whereas its imaginary parts at the bound-state masses are related to the decay 
widths. 



B. Higher n-point functions 

Dyson-Schwinger equations for higher n-point functions may be derived in a way that 
is similar to the case of the two-point function. Before showing them we need some more 
graphical rules (see Fig. 10), 
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n( P ) 



n n.f. 



T ( Pl +p 2 ,pj ,p 2 ) 




,n.f. 



R (Pi ,p 2 ,p 3 ,p 4 ) 



M (n) ( Pl ... p n ) 




Fig. 10 

where M^ n \ of course, should have n external (amputated) boson lines. 
For the three-point function e.g. we find the Dyson-Schwinger equation (in momentum 
space) 



M^( Pl +p 2 , Pl ,p 2 ) = (2^) 3 [m(P) m + m 2 (SP) C J Pl +p 2 ) 



+ m 2 (SP)Jp 1 )+m 2 (SP)Jp 2 )+m 3 (PPP)Jp 1 +p 2 ,p 1 ,p 2 )} 



(119) 



where m(P) m and m 2 (SP) m (p) are analogous to (107), (108) whereas the last term is given 
by 



m 3 (PPP) m (p 1 +p 2 ,p 1 ,p 2 ) = PiPjPkGn'Gjj'Gkk'Ti'fk'iPi +P2,Pi,P2) 



(120) 



and Tijk is just the exact three-point function for general chiral indices and external couplings 
equal to 1, see Fig. 11. 





+ 



Fig. 11 

The essential point is that again, may be reexpressed entirely in terms of non- 

factorizable n-point functions, namely 

M (3 \ Pl +p 2 , Pl ,p 2 ) = {2^fP i P j P k Yl ii ,{p 1 +p 2 )U jf { Pl )U kk ,(p 2 )- 



■ {pi'j'k' + Gi'lGj'mGk'nT^ n (pi + p 2 ,Pl,p 2 

)) (121) 
or, graphically (see Fig. 12; we ignore an overall factor {^.^fn) 3 on the r.h.s. of Fig. 12), 
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Fig. 12 

where the non-factorizable three-point function T n . f . is given by Fig. 13. 





+ 




+ 



Fig. 13 

The actual validity of (121), Fig. 12, has to be checked by a closer inspection of the 
Feynman graphs (it is just tedious combinatorics). 

We find that the non-factorizable n-point functions in our theory play a role analogous 
to the 1PI Green functions in other theories. 

The four-point function may be treated along similar lines. Again, the Dyson- 

Schwinger equation allows to express in terms of (P . . .) c m and (S . . ) c m n-point functions 
(which we show in coordinate space this time), 

M yi M y2 M y3 M y4 (HyiMy2My3My4)) c = 



16tt z 



m(S(y 1 ))8{y 1 - y 2 )S(y 1 - y 3 )S(yi - 2/4) + 



m 2 5( y y 1 -y 2 )5( y y 3 -y 4 )(S( y y 1 )S( y y 3 )) c + perm. + 



m 2 S(y 1 - y 2 )S(y 1 - y 3 )(P(y 1 )P(y 4 )) c + perm. + 



m 3 5(y 1 - y2)(S(y l )P(y 3 )P(y i )) c + perm. + 



m 4 (P( yi )P(y 2 )P(y 3 )P(y 4 )Y 



(122) 



Further, may be reexpressed in terms of non-factorizable n-point functions and 

reads 



M^( Pl , ...,p 4 ) = (47r) 2 P i P j P k P l U ul ip 1 )U jjl ip 2 )U kk/ {p 3 )U ll/ {p 4 ) \g vrk , v 



+Qi 'mQj'm'Qk'nQl'ri' 
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+ {Gi'j'm(ttmm'(Pl + P2) ~ $mm>)5m>k>l' + perm, J 
+ (ft'j' m n rom '(pi +p 2 )T^ nn ,(p 1 +P2,P3,P4)Gnk>Gn>l> + Perm. ) 

+ {pi'mGj'rn' T mm'n(Pl + P2, Pi, P2)Gnn'^n'r(Pl + P2)T^ S {P2, + Pi, P3, P^Qr'k'QsV + Perm. ) 

(123) 

where momentum conservation requires p\ + p 2 = Pz + Pa- Graphically, this identity may be 
depicted like in Fig. 14 (again we suppress an overall factor (2 v /7r) 4 in Fig. 14). 




^^^p vw H3~ w< ^^ + perm + 




+ perm 



The permutations in Fig. 14 contain all attachments of the external H{pi) lines that are 
topologically distinct (i.e. 3, 6 and 3 permutations). 

Observe that in each of the third type diagrams of Fig. 14 the lowest order diagram 
has to be subtracted in order to avoid an overcounting (this is so because IT(p) contains the 
lowest order, ^n(p) = Q + o(gl)). 
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X. BOUND-STATE MASSES AND DECAY WIDTHS 



A. General bound-state structure 



We claimed at the end of Subsection 9.1 that we could infer all the boson states and 
decay widths from the two-point function H(p), (117), which we want to discuss now. First, 
observe that the IT(p) propagator also occurs in higher bosonic n-point functions. E.g. for 
MW (Fig. 14), when one takes the third type of diagrams and inserts the lowest order 
(II ~ 1) for the four external H(pi) lines, there remains precisely an internal H(pi + P2) 
propagator (times Q). Therefore it is not a surprize that we can provide information on 
higher bosonic states, too, from Il(p). In fact, most of the information may be inferred 
from the denominator N(p), (118), of n(p). The zeros of the real part of N(p) will give all 
the bound-state masses of the theory - at least the leading order contribution - and the 
imaginary parts will give the corresponding decay widths ( 

This denominator N(p) reads, in lowest order 



62.6 



N(p) = 1 - a(g e ,p) - a{g* e ,p) 



(124) 



where, again in lowest order 

<x(ge,p) = 9eE + (p) , P(ge,p) = g$E-(p) , 
and the E±(p) are the exponentials of bosonic propagators, 



E ± (p) = J2(±l) n d n (p) , d n {p):-- 



(47r) n 



m—e 
2 



16 



(125) 



n=l 



nl 



Dlip). 



(126) 



The d n (p) are just n-boson blobs (see Fig. 15 for d 2 , d 3 ) 





Fig. 15 



and have the following properties: at s = —p 2 = (n/x) 2 , d n (p) has a singularity (real 
particle production threshold), and above this threshold it has an imaginary part. Therefore, 
slightly below the threshold (n/i) 2 , d n (p) is large enough to balance the coupling constant 
and make the real part of N(p), (124), vanish, 



rriE cos 6 d n (p) ~ 1 + o{m) 



(127) 



and, therefore, causes an n-boson bound state. At the position of the two-boson bound state, 
s = Mf = 4u 2 - A 2 , N(p) has no imaginary part and, therefore, the two-boson bound state 
is stable. At the three-boson bound-state mass M3, ^(s = Mf) has an imaginary part and, 
therefore, a decay into two Schwinger bosons (with mass u) is possible. For higher n-boson 
bound states the functions <i 2 , . . . , <i n _i have imaginary parts at M 2 , therefore decays into 
2, . . . , n — 1 Schwinger bosons are possible. 
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So far this is a lowest order reasoning, but we will find that we have to take into account 
some higher order effects, too, in order to obtain the physical spectrum of the theory. 

One higher-order effect may be understood easily. Remember that the a(gg,p) of (124) 
stems from the non-factorizable two-point function A, (114), (in fact, a is the ++ component 
of QA). To get more insight, let us rewrite A in terms of internal bosons, see Fig. 16. 

x^<+0+ ... +>A<+ ... ... 

Fig. 16 

We find that the one-boson line propagating from the initial to the final external vertex 
acquires no correction, because such corrections would be factorizable and are, therefore, 
excluded from a. On the other hand, all the higher d n (p), n > 2, do get corrections. 

This means that for the computation of the lowest pole mass (the Schwinger mass) 
one needs the bare Schwinger mass as an input, and the renormalized Schwinger mass is 
provided by the computation. For the higher bound states, on the other hand, one needs the 
renormalized Schwinger mass as an input in order to compute the bound-state mass poles. 
The reason is that the mass corrections for the bosons just shift the position of the threshold 
singularity and are therefore important in lowest order. There are other corrections present, 
too (internal boson interactions), however, they are unimportant in lowest order. 

This result is very plausible physically: the higher bound states should consist of physical 
Schwinger bosons with their physical masses \i (not the bare masses fi ). 

Therefore we redefine the d n (p) (n > 2) for the rest of the paper to be 

dn( P ) ■= { —rDji(p). (128) 

So we found, up to now, bound states composed of an arbitrary number of Schwinger 
bosons, where \i and M 2 are stable, and the higher bound states are unstable. 

However, this can not yet be the whole story. To understand why, look at the lowest 
order contribution to the three-point function, Fig. 12, with one incoming U(pi), one vertex 
and two outgoing U(p 2 ), n(p 3 ). Suppose the incoming H(pi) is at the mass — pi = M% of a 
sufficiently heavy unstable bound state. For a decay into stable final particles all the stable 
mass poles of IT(p 2 ), n(p 3 ) are possible. But by our above arguments the mass pole of the 
stable M 2 particle is present in IT(pj) as well as the mass pole of the Schwinger boson /i. 
Therefore, Fig. 12 describes decays into M 2 particles as well as fi particles. On the other 
hand, we did not find imaginary parts (up to now) in N(p) that describe decays into some 
M 2 , so obviously something is missing. 

The M 2 bound state itself was found by a resummation, therefore it is a reasonable idea 
to use the higher order contributions to a, (3 for a further resummation. a and (3 are just 
components of the non-factorizable propagator A{p), (114), so let us investigate it more 
closely. 

By a partial resummation we may find the following contribution to A(p), 

/cpq 
j^-^5 ijk A j f(q)Q j nU ll/ (q)Ai' k '(q)A km (q - p)k'k>m- (129) 
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This is just a blob where A(q — p) runs along one line, the other terms run along the other 
line. We want to discuss the [i-M 2 contribution, therefore we substitute A(q — p) by its 
lowest order, one-boson part, 

A(q - p) ~ AnD^q - p) ^ _^ j (130) 

and the two further »4(g) by their lowest order contribution 

The resummation that we need in (129) is taken into account by n(g). With these restrictions 
H{p) corresponds to the graph of Fig. 17. 




Observe that all internal bosons may be renormalized to their physical masses //, because 
this does not spoil non-factorizability in Fig. 17. The two factors A(q) in (129) are necessary 
in order to avoid an overcounting, but they cannot influence the presence of higher poles in 
Fig. 17. 

Now suppose that H(p) is at the M 2 + /i-threshold, s = —p 2 = (M 2 + fi) 2 . Then 
Dfi{q — p) is at its /i-singularity and 11(g) at its M 2 - singularity, and Fig. 17 corresponds (up 
to a normalization) to a ji-M 2 two-boson loop, i.e. Fig. 17 may effectively be substituted 
by Fig. 18, 



Fig. 18 



where the double line represents the two-boson bound-state propagator. 

Therefore, H(p) is singular at — p 2 = (M 2 + /i) 2 , and has a large real part slightly 
below and a large imaginary part slightly above this threshold. As a consequence, when the 
contribution of H(p) to a(p) is taken into account in the denominator N(p), (118,124), it 
will give rise to a further /i-M 2 bound state slightly below s = (M 2 + jj) 2 . Further it will 
open the /i-M 2 decay channel at s = (M 2 + /i) 2 . 

Now suppose we put 11(g) in (129) on a higher (unstable) bound-state mass M n , n > 2. 
Then in the denominator N(q) of 11(g) the real part again vanishes, but there remains an 
imaginary part. Therefore, H(p) is finite and imaginary at s = —p 2 = (M n + /i) 2 and cannot 
give rise to a fi-M n bound-state formation. 

Further, because there is no threshold singularity at s = (M n + /i) 2 , this means that no 
new decay channel opens at that point (i.e. the imaginary part of H{p) varies smoothly 
around s ~ (M n + /i 2 )), which simply means that the unstable higher n-boson bound states 
are no possible final states (of course, they are possible as intermediate resonances). 
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We could substitute the one-boson line in Fig. 17 by another AQUA line and would find 
that this graph behaves like a M 2 -M 2 blob near s = (M 2 + M 2 ) 2 , and we could allow for even 
more AQUA lines. The physical picture that evolves from these considerations is like follows: 
in addition to the unstable n-boson bound states there exist further (unstable) bound states 
that are composed of Schwinger bosons \x and (stable) two-boson bound states M 2 . Further, 
the unstable bound states may decay into all combinations of p and M 2 particles that are 
possible kinematically. The imaginary parts of the corresponding n-particle blobs (where 
particle means p or M 2 ) are large near their thresholds, therefore there is a kinematical 
tendency to rise the decay probabilities for decays with small kinetic energy. This is not so 
surprizing, because in 1 + 1 dimensions the phase space "volume" does not grow with kinetic 
energy. 

Further we want to emphasize the following point for later convenience. We found, by 
a further resummation, contributions to N(p) that may be substituted by a n-particle blob 
that contains p and M 2 bosons (or only M 2 ), near their respective n-particle thresholds. 
This is true for the real parts. The imaginary parts are given precisely by the threshold 
singularities which stem solely from the n-particle blobs. Therefore, the imaginary parts of 
these resummed contributions may be substituted by the corresponding imaginary parts of 
the n-particle blobs for arbitrary values of p. 

Finally let us briefly comment on the special case = 0. Here parity is conserved 
and we have to investigate scalar (SQU(p)S) and pseudoscalar (PQYl(p)P) propagators 
separatly. We find a partial cancellation between numerator and denominator in (117) 
(a(ge=o,p) = a(g% =0 ,p), etc.) 

SQU(p)S "' E 



1 - a(ge =0 ,p) - (3(g e=0 ,p) 



PQK(P)P = z ? r ^- W( r (132) 

1 - a{g e=0 ,p) +p{ge=o,P) 

and, therefore, using the lowest order expression (125), we find that the odd (even) n-boson 
blobs d n (p) are cancelled in SQU(p)S (PQTl(p)P), so that only even (odd) mass poles remain. 
Of course, these parity considerations may be generalized easily to the mixed bound states 
(where each p is odd and each M 2 is even). 

B. n-boson bound-state masses 

The mass pole equation of the n-boson bound state in lowest order is given by (see (127)) 

f n {p) := 1 - mS cos dd n {p) = (133) 
For the Schwinger boson (n = 1) this equation reads 



-47T 

P 2 + Vo 

with the solution 



l = mEcos#— 2 (134) 
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— p 2 =: /if = /il + AnmH cos 6 —: + Ai 



(135) 



(we will denote the lowest order corrections to all the pole masses M 2 by A n ). 

When we want to recover the second order result for the Schwinger mass of Section 8 we 
have to include all terms of N(p), (118). The mass pole equation up to second order reads 



-47T 
P 2 + 



(9i + 91 + 92+ 9* 2 ) + {91 + gl)E^\p)- 



9i9i 



(E?( P ) + 



V KJ P 2 + Vo 



-An 



v 2 -vl' 



(136) 



where g± (#2) are the first (second) order contributions to gg (see (89)) and E±\p) are the 
exponentials E± (p) without the one-boson term, see (52). In (136), the g^E^ (p)) 2 parts 
are, in fact, o(m 3 ), and may be omitted. The solution to (136) is 



2 22 
H 2 := -p = fi Q 



1 + 4vr^ cos 9 + 2vr^^((£ + + ^(1)) cos 29 + E_ - E ( }\l))} (137) 



and, indeed, coincides with (96) of Section 8 that was obtained by a direct perturbative 
calculation. 

For the two-boson bound state mass in lowest order we have to solve 



1 



l=y(gi+gl)16n 2 (Dl)(p) 



(138) 



where now fi is the physical Schwinger mass (137) including fermion mass corrections. (D 2 ) (p) 
is just the two-boson blob of Fig. 15 and may be evaluated by standard methods: 



(Dl)(p) 



d 2 q -1 



-1 



(2vr) 2 q 2 + /i 2 (q-p) 2 + /1 2 



r d 2 q r 1 
J Jo 



dx 



(2tt) 2 Jo [q 2 + 2pq(x - 1) + p 2 (l - x) + /1 2 ]' 



1 r 1 



dx 



1 /-i dy 



An Jo p 2 x(l-x)+/i 2 ix(-p 2 ) Jo yt + ^-l) 



1 1 1 

■ arctan 



ir(-p 2 ) R(p 2 ) 



R{p 2 ) ' 



(139) 



R(P 2 ) 



! Afi 2 
—p 2 



- 1 



(140) 



where we used the fact that, for the bound state, — p 2 has to be beyond the threshold, 
-p 2 < A/i 2 . 
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Now we simply have to insert this result into (138) in order to get the mass pole: 

SvrmEcos^ 11 

For small fermion mass m R(p 2 ) is very small, too. Therefore, the leading order result 
will stem from a matching between these two factors, where we may set — 4^2 and 
arctan -^^y = § • Doing so we get 



9 . 7r 2 mS cos 6 ,„ 
^(P 2 ) = 3 (142) 



or 



l .-, 7r 4 m 2 E 2 cos 2 6* 
M l := ^ \ { , ^ r e f - V(l ^ ) =: V - A 2 , (143) 



9 ^ 
777 777 

ilff = 4/i 2 (l — 7.83-y cos 2 6 + o(— y)) (144) 
^0 ^0 

which is of second order in m. Again, this result coincides with the one from a direct 
perturbative calculation ( f34|). 

For the computation of the three-boson bound-state mass M3 we need the three-boson 
propagator d 3 and find, in lowest order (see (133)) 

1 = ^mS cos 9 ■ 647r 3 L>3 (p) (145) 
or, after a rescaling p — > - to dimensionless momenta 



647r 3 mS 
6 /i : 



1 = — — CO sQDl(p). (146) 



D^(p) is given by the three-boson loop of Fig. 15 (in the sequel we introduce positive squared 
momentum s = —p 2 > 0) 



d 2 q x d 2 q 2 



rl rx r 

-2 dx dy 
Jo Jo J 



(2tt) 4 (p + g 1+g2 )2 + l^ + lg2 + l 

d 2 q 1 d 2 q 2 1 



q\ + 1 + {q\ - q\)x + ((p + gi + g 2 ) 2 - g|)y 



(47r) 2 7o s(a;|/ — x 2 y — y 2 + xy 2 ) — x + x 2 — xy + y 2 
dx [% dz 



]tt 2 (1 - s(l - x)) Jo z 2 + T 2 (s,x) 
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arctan — — -, (147) 



o 8tt 2 (1 - s(l -x))T(s,x) 2T(s,x) 

T 2 (s v _ x*-sx*(l-x)+4x(l-x) 

lS ' ] ~ 4(s(l-a;)-l) ' 1 j 

where, as usual, we introduced Feynman parameter integrals and performed the momentum 
integrations. Further, the first Feynman parameter integral could be done analytically. The 
numerator of T 2 has a double zero at s = 9: 

9x(x- 2 -) 2 . (149) 

This double zero is in the integration range of x and is precisely the threshold singularity. 
Setting 

s = 9-A 3 (150) 

in the numerator of T 2 in the factor ^, and s = 9 everywhere else, where it is safe, one 
arrives at: 



1 -i j arctan -^^3 

127T 2 7o ^| 9rr _ 8 | ^ _ 1)2 X + ^ X 2(! _ 

The mass-pole equation reads 

647T 3 

1 = — mS cos 0J(A 3 ) (152) 

and must be evaluated numerically. It gives rise to an extremely tiny mass correction A 3 . 
For sufficiently small m it is very well saturated by 

263 

A 3 (mS cos 9) ~ 6.993 exp( -) (153) 

mL cos & 

and is therefore smaller than polynomial in m. (I checked the numerical formula (153) 
for 30 < msLe < 1000 ' corresponding to 10~ 2 < A 3 < lO" 100 , but I am convinced that 
it remains true for even larger mS ^ osg ; however, there the numerical integration is quite 
difficult because of the pole in (151).) 

We conclude that the three-boson bound state mass is nearly entirely given by three 
times the Schwinger boson mass (we change back to dimensionfull quantities now), 

a 2 

M\ = 9/i 2 - A 3 , A 3 = 6.993^ 2 exp(-0.263— f -) (154) 

mS cos 

or, differently stated, that the binding of three bosons is extremely weak. 

Therefore it holds that M 3 > \x + M 2 , and, consequently, a decay of M 3 into fx + M 2 is 
possible. This has the consequence that the three-boson bound state is unstable even for 
9 = 0. 
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C. A mixed bound-state mass 



For our further discussion we will need the residues of the propagator Tl(p) at its various 
mass poles. The n-boson mass poles were, in leading order, the zeros of the functions / n (p), 
see (153). Therefore, we need the first Taylor coefficient c n of each f n (p) at its mass pole 
M 2 n = (n/i) 2 - A n (n > 2) or M 2 = $ + A x , 

fn(s) — c n(s — M 2 ) , C n = ^j-f n {s)\ s=M 2 (155) 

where s = —p 2 . The c n may be easily obtained from our mass computations. From (134) 
and (141) we find for c\ and c 2 

Cl = AnmLosO = A^ (156) 



° 2 8vr 4 (mScos^) 2 2A 2 (157) 

For the computation of C3 we observe that because of formulae (150), (153) mS cos 9ds(s) 
may be written, in the vicinity of s = Mf , like 

v qa ( \ mScos #i 6.993/i 2 

mE cos 9 ds(s) ~ In — . 158 

3V ; 0.263 V - s y ' 

Therefore, we find the Taylor coefficient 

m£cos# 

C3 = 0263^- (159) 

Further we need, for the computation of the lowest mixed bound-state mass M^i (which 
is composed of one M 2 and one ji = Mi), the residues of the propagator H(p) at the two 
lowest mass poles. The denominator N(p), (118), is given by (155), and for the numerator 
of n(p), (117), we use a(g g , s ~ M 2 ) ~ geE + (s ~ M 2 ) ~ ^x^, etc., which holds near the 

pole, see (124) - (127), and find (here n = 1,2) 

1 J (^ + ^)c„( s -m 2 ) U-ir^ sg ; libUj 

For the computation of the p,-M 2 bound state we need, in addition, the matrix A, (113), at 
the n-boson mass poles (here 1, because there A itself has a pole like (160)), 

Now we are prepared for the computation of the M 1;1 mixed bound-state mass. In Subsection 
10.1 we claimed that we could substitute the resummed contribution H{p) to N(p), see (129) 
and Fig. 17, near its threshold by a two-particle blob consisting of one u and one M 2 (times a 
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not yet specified factor), see Fig. 18. This we achieve by inserting expression (160) for H(p) 
near its two poles /z, M 2 , into (129), and by using (161) for the remaining A(—q 2 ~ Mf). 
Althogether we find for H(p) near its threshold 

H iV {- P 2 ~ ( m 2 + /i) 2 ) ~ / * yfc _L_^ (]]) ( 9 ° °A r—^r-2 Tm - 

J (2tt) 2 g e + g* \ 1 1 / jy V #e (<?e + 9*e)c 2 (-q 2 - Mf) 
f 5 , g* g \ 1 ( 1 1\ 1 f 1 -1 ^ 

Oi'k'm 



£e 9o ) w 9o + 9o \ l 1 / w + ^)ci (-(<?- p) 2 - /i 2 ) 



/ A i ; 



^0 O^" 1 l 1 ),/^'/^^ (162) 
The contribution of H u i to ct(^) + cx(gg) in the denominator N(p), (124), is 

9eH++{p) + g* e H—{p) =: (g e + g* e )d 1 , 1 (p) = 

• d 2 q 8vr 4 m£cos# 4tt 



(9o + 9*o)Jj^ 



(2tt) V(g 2 + Mf) (p-q) 2 + fx 2 
327r 5 m 2 S 2 cos 2 6 



arctan 



27r / u 2 w(s,M|,/i 2 ) 
2s 



(vr+ 



w(s, Mf, v 2 ) - H J W) {s + fi 2 - M 2 )(s -f, 2 + Mf) 



) (163) 



w(x, y, z) := {—x 2 — y 2 — z 2 + 2xy + 2xz + 2yz) 2 (164) 
where s = —p 2 ~ (/i + M 2 ) 2 . The /i-M 2 bound-state mass fulfills the equation 

l = {9o + 9e)diM (165) 
with the solution in leading order (here M^i denotes the /i-M 2 bound-state mass) 

,,2 / > o . . 327r 10 (mScos^) 4 , , 

M 2 hl = (/i + M 2 ) 2 - Ax,! , A 1)1 = ^ ^ (166) 

which is valid for sufficiently small A^. M^i was computed from a two-boson blob (Fig. 
18), like M 2 , therefore the Taylor coefficient of (s — M 2 X ) is analogous to c 2 , equ. (157), 

1 nfi? x 

01,1 2Ai,i 647r 10 (mEcos^) 4 ' 1 j 

Further, the above equ. (162) shows that the /z-M 2 blob enters into the functions a, 

/3 of „4, (113), like any other odd n-boson blob d n {p). 

In principle, even higher mixed bound-state masses could be computed along similar 
lines, but we want to change now to the computation of the decay widths of the lowest 
unstable bound states. 
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D. Decay width computations 



In order to find the decay widths of some bound states, we have to examine the imaginary 
parts of the denominator N(p) in the vicinity of the corresponding mass poles. Using the 
first order approximation (125) for a, [3 we find for N{p) ( e.g. in the vicinity of s ~ Mf for 
definiteness) 

m 2 T? 

N{p) ~ 1 - mE cos 9E + (p) + -^—{E 2 + {p) - E 2 _(p)) 

= 1 - mE cos 9 (dx(p) + d 2 (p) + d hl (p) + d 3 (p) + ...) + 

m 2 T 1 2 ( y d 1 (p){d 2 (p) + d 4 (p) + ...) + d ltl (p)(d 2 (p) + d 4 (p) + ...) + 

d 3 (p)(d 2 (p) + d^p) + ...) + ...) (168) 



where we included the [i-M 2 blob d^i, as discussed above, because we need it for the subse- 
quent discussion (we ignore, for the moment, higher M 2 blobs that are, in principle, present). 
Near s = M 2 the real part of (168) is given by c 3 (s — Mf ) and we find 

N(s ~ Mf) ~ c 3 (s - Mf) - imll cos 0(lmd 2 (s ~ M 3 2 ) + Imd 1;1 (s ~ M 3 2 )) + 



im^E 2 d 3 (s ~ M3)Imd 2 (s ~ M 3 ) + o(m' 
c 3 (s - M 3 2 ) - imT,(cos9 ^-)Imd 2 (M 3 2 ) - imT, cos 9 lmd 1A (Ml) + o(m 2 ) 

COS f/ 

(169) 



where we used d 3 (Ml) ~ mScosg , see (133). 

This computation may be generalized and tells us that parity forbidden imaginary parts 
(decay channels) acquire a factor (cos^ — ^^), whereas parity allowed imaginary parts have 
the usual cos 9 factor. 

[Remark There seems to be something wrong with the sign of the parity forbidden 
imaginary part (the d 2 term). Actually the sign is o.k. and the problem is a remnant of the 
Euclidean conventions that are implicit in the whole computation (see Section 1). In these 
conventions 9 is imaginary and therefore (cos# — > 0. Of course, this is not a reasonable 
convention for a final result. When performing the whole computation in Minkowski space 
and for real 9, roughly speaking, the roles of E + and E- are exchanged in (168). This gives 
an additional relative sign between parity even and odd n-boson propagators and, therefore, 
changes the factor of d 2 to (^^ — cos#), which is > for real 9. We will keep this remark 
in mind and express the final results in Minkowski space and for real 9] 

Now we may find the M 3 decay widths by comparing the inverse of (169) to the general 
formula 

n< \ const. 

G(P) ~ (170) 



where T is the decay width. We find (for real 9) 

1 const. 



N(s ~ Mf) s-M 2 - - cos0)rmd 2 (M 3 2 ) + cos 9\m d 1;1 (M|) 



(171) 
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Next we need the imaginary parts Im^, Imd^i. Both of them stem from a two-boson blob, 
so let us write down the general result (which is standard) (s = — p 2 ) 



lm(D Ml D M2 )(s) = Im J 



d 2 q -1 -1 

(2n) 2 q 2 + M 2 (p-g) 2 + M| 
1 



2 W ( S ,M 2 ,M 2 )' 



(172) 



w(x,y,z) = (x 2 + y 2 + z 2 - 2xy - 2xz - 2yz) 2 . (173) 

2 

Therefore we can write for (171) (where the normalization factors of d,2 and d lt i are ^ and 
Cic 2 , respectively) 



const. 



(174) 



o _ A/f 2 _ jmZ( 1 pnn m 4tt 2 _ ■ (mS cos 6>) 2 l&nj 

and therefore, by using the approximations 

w(M 2 , /i 2 , /x 2 ) ~ w(9fi 2 , fi 2 , fi 2 ) = zVE/i 2 (175) 

w(M 2 , M 2 , fi 2 ) ~ w(9/i 2 , M 2 2 , /i 2 ) = 2a/3a* V^3 + °(™ 2 )> ( 176 ) 
we find the following results: 

Tm 3 ^ = 0.263^^(^- - 1) 
9V5/i cos 2 6 

~ 3.608//(— ^ - 1) exp(-0.929— L) (177) 
cos 2 9 ' mcosfc* 



and 



4tt 3 A 

^M 3 -*fi+M 2 — 0.263 



3v^3/i 

~ 43.9uexp(-0.929 — ^— ) (178) 
mcosO 

where we inserted the numerical value E = ^ = 0.283/x. 

The ratio of the two partial decay widths does not depend on the approximations that 
were used for the M 3 computation, 

r _^ 1 

M 3 ^2/i _ C; : j« (179) 



rM 3 ^+Af 2 V157T 

Analogously we may compute the decay width of the mixed bound state M^i, starting from 
iV(s ~ M^) ~ ci,i(s - MjJ - im£(— !— - cos ^)Imrf 2 (M 1 2 1 ) (180) 

COS (7 
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which leads to the decay width 



2 8 7r 12 (m£cos#) 5 / 1 , ^mcos^,, 1 . . . 

r — 2 " = ( ^ - l) K 21340 " ( — » ( ^ - 11 (181) 

for the decay M^i — > 2//. This decay is parity forbidden, and therefore Ml^ is stable for 
6 = 0. ' 

In principle, we could have computed the above decay widths by another method, too, 
namely by the use of the resummed three-point function of Fig. 12. Choosing the first 
graph on the r.h.s. of Fig. 12 (consisting of three exact propagators and one pure vertex), 
we could precisely rederive our results (177), (178), and (181). 



XI. SCATTERING 

A. Two-dimensional kinematics 

For a discussion of scattering processes we need some basic facts about two-dimensional 
kinematics. We will restrict our discussion to elastic scattering. Suppose we have two 
incoming particles with masses M 1; M 2 and momenta p±, p 2 , and two outgoing particles, 
again with masses Mi, M 2 , and with momenta p 3 , p±. Momentum conservation requires 

p ■=p 1 +p 2 =p 3 +p 4 (182) 
and all momenta are Minkowskian in the sequel. In the center of mass system we may write 



(Jk 2 + MA (Jk 2 + M. 

Pl= [ k 



, Jk 2 + M 2 \ A/A; 2 + Mf\ 

P3 ={ ±k ) ' P4= { T k ) (183) 

where in p 3 , p 4 the first sign is for transmission, the second sign is for reflexion. For the 
kinematical variables we find for transmission 



s = (pi + P2) 2 = Ik 1 + Ml + M 2 + <lJ{k 2 + M 2 )(k 2 + M 



t T = ( Pl - p 4 ) 2 = -2k 2 + Mj + M 2 - 2\J{k 2 + M 2 ){k 2 + Mf) 

u T = (pi ~ Ps) 2 = (184) 



and for reflexion 



s = (pi + P2) 2 = 2k 2 + Ml + Ml + 1\J (k 2 + M 2 )(k 2 + Mf) 

*fl = (Pi - P 4 ) 2 = 2A; 2 + Ml + M| - 2 V /(F + M 1 2 )(A; 2 + M 2 2 ) 

«i? = (pi-p 3 ) 2 = -4A; 2 (185) 
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When the two masses are equal, the two particles are identical in our theory and the discrim- 
ination between transmission and reflexion does not make sense. The kinematical variables 
turn into 

s = (Pi + Vif = 4(A; 2 + M 2 ) 
t = (Pi - P±? = -4A; 2 

u=(pi-p 3 ) 2 = (186) 
The elastic scattering cross section of two particles is given by 

,x C syni \M(s)\ 2 

where u> is defined in (173), .M is the transition matrix element and C sym is a symmetry 
factor that takes into account identical particles in the final state (C sym = j 1 , for n x 
particles M x and n 2 particles M 2 in the final state). As it stands, expression (187) holds 
provided that the initial and final particle propagators are normalized in the usual fashion 
(~ s_mO - Otherwise, (187) is multiplied by the normalization factors (the residues of the 
propagators) . 



B. Scattering processes 

Finally we are prepared for a discussion of scattering. Let us focus for the moment on the 
lowest order graph of Fig. 14 for the four-point function (123). It consists of four external 
exact propagators H(pi) and a simple vertex as the lowest order transition matrix element. 
The n(pj) contain two stable-particle mass poles, \x and M 2 , therefore this graph describes 
/i and M 2 scattering (this remains true for higher order contributions; as a consequence, the 
same transition matrix elements contribute to /i and M 2 scattering processes, and they may 
only differ by some kinematical and normalization factors). 

Let us consider elastic scattering of two Schwinger bosons for definiteness. Then each 
external H(pi) propagator is odd and contributes to the graph like (s, = — pf) 

n #(s , = ,V, = ^i#( 1 1 ) (188) 

s i f 1 \ V j 

Here we face the problem that the first graph of Fig. 14 is already of fifth order, because 
each propagator II(sj) has an external vertex. We just omit these external vertices (i.e. we 
omit the factor (ge + gg) in (188) for each propagator), because we want to discuss first order 
scattering. Doing so, we find for this graph 

4 47T 4 47T 

Ph p h P 33 P 3^hhkiGkik^iziik2 II _ 2 = (9° + 9e) II _ 2 ( 189 ) 

i=l S i A* i=l S i f 1 

i.e. each /i propagator has a residue 4n. In order to obtain the transition matrix element 
one has to amputate the external boson propagators in the usual LSZ fashion. When the 
propagators are normalized by r\ = Att, the bosons themselves are normalized by v^4vr, which 
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has to be divided out for each amputation. This leaves a factor V^r for each external boson 
in the transition matrix element. However, the squared transition matrix element enters the 
scattering cross section, therefore the net effect on the cross section is a multiplication by 
the corresponding propagator residue for each external line. 

Therefore, we find for the lowest order boson-boson elastic scattering 

. . „ |(mS cos 6) 2 

wW'H ^^j (wo) 

where we have for the propagator residues (they may be inferred from (160), r\ = c .^ +g *^ ) 

87r 4 m£cos# 

r 1 = 47r , r 2 = . 191 

(190), of course, coincides with a naive computation using the first order bosonic four-point 
function ($(xi) . . . <&(x4)) c m (the latter may be inferred immediately from (48)). Observe 
that (190) is singular at the real particle production threshold s = 4/i 2 (w(4/x 2 , fj, 2 , /i 2 ) = 0). 

In a next step we want to consider the second order contribution of Fig. 14 (the third type 
graphs). There are three graphs of this type, namely s, t and u channel, but we will consider 
only the s channel (annihilation channel) for the moment. In this diagram the lowest order 
graph must be subtracted in order to avoid overcounting (see Fig. 14), therefore the graph 
of Fig. 19 




Fig. 19 



contains the lowest order graph and the second order s channel contribution. 
Actually we will allow for arbitrary final states in the sequel, // + //—>•/, because this 
enables us to use the optical theorem, which may be written for the current problem like 



ff "-/ (s) = ^Mm lmM ^ [s) (192) 



where M a b^ab is the forward elastic scattering amplitude. w(s, M 2 , M 2 ) is an initial state 
velocity factor; the final state factors must be produced by M. a b^ a b, as we will find in the 
sequel. 

Specifically we choose a — b — /i, and, therefore, both vertices of M. are contracted by 
scalars S (we use matrix notation) 

M 2 ^(s) = S T QU(s)S (193) 

Before starting the computations, we want to make some comments. First, as is obvious 
from Fig. 19 and our discussion, in (192) all combinations of ri\ji and n 2 M 2 are allowed as 
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final states. Consequently, they must exist as intermediate states in M. 2il -^ 2il , too, in order 
to saturate the optical theorem (192). Therefore, we are forced to include the M 2 particle 
into the two-point function IT(p), as we did in the previous section, in order to maintain 
unitarity. 

Secondly, in finite order perturbation theory the optical theorem relates graphs of differ- 
ent order. However, we use a resummed perturbation series in (192) and, therefore, we will 
find a relation that holds for the whole, resummed two-point function IT(s). 

In a first step we want to discuss the special case 9 = 0, because it is much easier and 
shows the relevant features without technical complications. For 9 = the amplitude (193) 
reads (see (132)) 

M e=0 (s) = = 

' 1- ?f(E + (s) + E_(s)) 

^ (194) 



1 - mE(d 2 (s) + d 2 fi(s) + d 4 (s) + . . .) 

where we inserted the lowest order (125) and expanded the exponentials E±(s) like in (126). 
Again, we include the M 2 particle (which is found by a further resummation) into E±, 
because this is absolutely necessary, as we have just argued. Actually d 2t0 describes the M 2 - 
M 2 blob, and in (194) only parity even contributions may occur. For the optical theorem 
(192) we need the imaginary part 

TmA// e=o /x m 2 £ 2 (Imrf 2 (g) + Imrf 2 ,o(s) + lmd 4 (s) + . . j 

1 [1 - mL{Red 2 {s) + . . .)] 2 + m 2 £ 2 (Imd 2 (s) + . . .) 2 1 ' 

We find the following physical picture: at s = 4/i 2 the elastic scattering threshold (/ = 2/x) 
opens, at s = 4M| the 2/j, — > 2M 2 threshold is added, at s = 16/z 2 the 2/x — > 4// threshold, etc. 
The d n (s) were defined as d n (s) = ^D™(s), therefore their imaginary parts are precisely the 
final state factors for the corresponding cross section, including the phase space integration 
(the cutting of the D"(s)), the propagator normalizations r\ = 4n, and the final state 
symmetry factors for n identical particles, C sym = For the multi-M 2 propagators d m ,o(s) 
(and, more generally, for d mjJl (s)) the first two points (propagators with their residues) are 
obvious, the third one (correct final state symmetry factor) may be checked by a closer 
inspection of the mass perturbation series. We show it for the lowest order contribution to 
the M 2 -M 2 propagator d 2;0 (s), where we depict in Fig. 20 this lowest order contribution and 
the perturbation expansion graph where it stems from. 



Fig. 20 

The second graph in Fig. 20 is a second order mass perturbation, therefore it contains a 

2 

factor ^j-. Further there exists precisely one diagram of this kind in the perturbation series, 
therefore the ^ factor remains in d 2>0 (s) as the required final space symmetry factor. Via 
some combinatorics this argument may be generalized to higher order contributions to the 
M 2 -M 2 loop d 2:0 (s) and to higher multi-M 2 loops. 
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For the total cross section (192) we get 



tot,0=o/ \ r 2 m 2 E 2 (Imd 2 (s) + Imd 2 ,o(s) + Imd 4 (s) + . . .) 

" } w(s,/i 2 ,/i 2 )([l-mS(Rerf 2 (s) + ...)] 2 + m 2 S 2 (Imrf 2 (s) + ...) 

^ = Im^T 5 ) etc ' (197) 

which we want to evaluate for some specific values of s. At the elastic scattering threshold 
s = 4/z 2 , lmd 2 (s) is singular and we find 

<i/ =0 (V) = 4. (198) 

Therefore, the singular behaviour of the lowest order cross section at s — 4/x 2 is cancelled 
by higher order contributions. This behaviour is, however, further changed by the t and u 
channel contributions. 

In an intermediate range, far from all thresholds and bound state masses, 4/z 2 < s < 4M|, 
cr tot is well described by the lowest order result (190), because there mEd n (s) is small 
compared to 1, 

toti , =0 ^ rlm 2 mmd 2 (s) _ fcjmg 

2 ^ / [ } ~ w(s, /i 2 , /i 2 ) ~ 2w 2 (s, /x 2 , /x 2 ) ' 1 

At the first bound-state mass, s = Mf < 4M|, a resonance occurs. There the real part 
contribution to the denominator of (196) vanishes by definition and we find 

tot 0=0/ 2 x_ r 2 m 2 £ 2 Imd 2 (M 2 ) 
<V*/ ( M 2,o) = w (M 2 ,/x 2 , / i 2 )m 2 S 2 (Imrf 2 (M 2 )) 2 = 4 (2 ° 0) 

and the resonance height does not depend on the coupling constant (of course, the width 
does). 

At the 2M 2 production threshold s = 4M| the scattering cross section goes down to zero 
(here d 2 $ is singular) 



totfi=Q 2 _ r\m 2 Y? Imrf 2 , (4M 2 ; 

2 ^ 1 2) ~ w(AM 2 ,ii 2 ,ii 2 ) m 2 S 2 (Im4,o(4M 2 )) 2 1 j 



In addition, at this point the 2/j, — > 2M 2 production channel opens. At the four-boson 
bound-state mass s = M\ we find the next resonance 

tot,e=0/ M 2, = r 2 (lmd 2 (M 2 ) + lmd 2fi (M 2 )) 

2 ^ f 1 4) w(M 2 ,fi 2 ,fi 2 )(lmd 2 (M 2 ) + Imd 2 , (M 2 )) 2 {U) 

Again, the resonance height does not depend on the coupling constant, and, in addition, 
here already two decay channels are open for the M 4 resonance. 

At the 2/j, — > 4/i real production threshold s = 16/i 2 , a tot again vanishes, and for even 
higher s the above pattern repeats. 
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Observe that, because <7 tot has a local maximum (resonance) at the bound-state masses, 
whereas it is zero at the real particle production thresholds, the resonance widths (decay 
widths) must be bounded by the binding energies. For the M 1;1 and M 3 decay widths this 
may be seen from the explicit results (177), (178) and (181). 

The t and u channel contributions do not change this pattern (they have no imaginary 
parts and are small for all t, u). 

Next let us turn to the ^ case. There parity forbidden transitions are possible, and 
therefore we will find M ltl and M 3 resonances, too. The forward scattering amplitude (193) 
reads 

M (s) = 9e + 9l-j9e9l{E + {s)-E_{s))^ 

^ 1 -(9o + 9*o)E + (s) + gog*e(El(s) + E 2 _(s)) " 



9e + 9* e - 4ffgffg(di(s) + d hl (s) + d 3 (s) + . . .) 

l-(ge + g* e )(di(s) + d 2 (s) + d hl (s) + ...)+ ±g e g*e[di{s){d 2 {s) + d 2>0 (s) + ...) + ...] 

Please observe the presence of only odd di in the numerator and of only odd x even di x dj 
in the second term of the denominator. Therefore, the ^gog% terms in the numerator and 
denominator do not contribute to parity allowed transitions, and the discussion of such 
parity allowed transitions is analogous to the 9 = case that we discussed above. 
Again, we want to discuss the scattering cross section 

<w< s > = ^k?) lmM ^ (s) (204) 

for some specific values of s. At s = 4fi 2 we find again 

tot Un 2, = r\ (gg + g* e ) 2 \md 2 {^ 2 ) 

2 ^ f{ } w(4^,n\^)l + (ge + g* e ¥(lmd 2 (4^))* { U ° j 

At the first parity forbidden resonance s = M\ x we find (Red^i = g9 + g * ) 

^tot (M 2 , rj go + g* e ~ 4»*Red 1 , 1 (M 2 1 ) 



u>(M 2 1; fi 2 , /i 2 ) -i(g e + g* e )lmd 2 (M 2 1 ) + 4^Red 1 , 1 (M 2 1 )Imd 2 (M 2 1 ) 

= 4 (206) 

and, therefore, the same resonance height as for the first parity allowed resonance in the 
9 = case (200). 

At the parity forbidden threshold s = (M 2 + /i) 2 , where lmd 1:1 ((M 2 + /i) 2 ) is singular, 
we find 

o:;:: . f i(M 2 ■ fif) 7 
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Im g e + g* e - Mg e g* e lmd hl ((M 2 + /j,) 2 ) = 

1 - i(g e + g* e )(\m.d 2 + Inic^i) - Ag e g* e \m.d 2 \m.d 1) i 

r\ (g e + g* e ) 2 (lmd 2 + lmd hl ) - 4g e g* 9 lmd hl (l - 4g e g* e lmd 2 lmd hl 

W ((M 2 + /U ) 2 ,/i 2 ,/i 2 ) " (l-4g e g* e lmd 2 lmd 1 , 1 ) 2 + (ge + g* e ) 2 (lmd 2 + lmd 1 , 1 ) 2 

r 2 Ag e g* e (lmd hl ) 2 lmd 2 , Ag e g* B x 2 rjlmd 2 ((M 2 + fif 



( *yey 6 y 

V A/, 4- nJ / 



W ((M 2 + / i) 2 ,/x 2 , y u 2 ) (^ + ^) 2 (Imd 1)1 ) 2 ^ge + gp w((M 2 + fi) 2 , fi 2 , fi 2 ) 

(207) 

where we performed the limit Imrfi i — > 00 and kept only the lowest order contribution in 
gg. Therefore, in contrast to the parity allowed case, the parity forbidden thresholds do not 
give zero in a tot . 

The reason for this behaviour may be easily understood. In the limit of 9 — > there 
should not remain any effect of resonances or thresholds in <7 tot for parity forbidden transi- 
tions, and <7 tot should be described by the lowest order result (190). 

Precisely this happens: Although the resonance height at M\ x remains unchanged for 
9 — > 0, (206), its width tends to zero, (181). This means that the resonance M^i still exists 
but is stable against M ltl — > 2/i decay for 9 — > 0. Actually the M ltl bound state is a stable 
particle at all for = 0. Further, at threshold s = (M 2 + /i) 2 , <r tot tends to the first order 
result (190) for 9^0, 

lim(i^) 2 = m 2 S 2 + o(m 3 ) (208) 

as it should hold. 

For even higher s, when both parity allowed and parity forbidden final states are possible, 
we again have the problem that the relative sign of the parity forbidden process is "wrong" 
due to our conventions (see the remark after equ. (169)). E.g. at the M 3 resonance we find 
from (203) 

tot 



g e + g* e - 4g e g* e Red 3 (M 2 



Im 



-i(g e + g* 9 )(Imd 2 (M$) + Imdi,i(M 3 2 )) + 4ig e g* e Red 3 (Ml)Imd 2 (Ml 

^96 9 „ 



+ 9o 



90 +9 2 



w(M 2 , v 2 , /i 2 ) {gg + g * _ ^l)l m d 2 (M 2 ) + (g e + g* )lmd 1A (M 2 ) 



mS(cos0- -±3) 

^ cos a 1 



w(MI,/j, 2 ,/j, 2 ) mS(cos0 - -Mrmd 2 (M 3 2 ) +mS cos 91md 1A (M$) 



cos 8 ' 

r? — — COS6 1 

I cos 9 



w(M 2 ,fi 2 ,fi 2 ) (^ - cos9)lmd 2 (M 2 ) + cos 9Imd 1A (M 2 ) 
r\ sin 2 #(sin 2 9lmd 2 (M 2 ) + cos 2 9\md ly i(M 2 )) 
w(M$, n 2 , /x 2 ) [sin 2 9lmd 2 (Ml) + cos 2 9lmd l , 1 (M 3 2 )] 2 



(209) 
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where the last two lines are for real 9 (the last line may be easily checked by a low order 
reasoning, in case somebody does not trust our imaginary 9 convention). 

Again, the resonance height (containing two partial decay channels) does not depend on 
the coupling constant. 

For even higher s the above pattern repeats. 

The last thing to be discussed is the contribution of the t and u channel diagrams. There 
the lowest order diagram must be subtracted (see Fig. 14), 

M' 2 ^(t) = S T g(U(t) - 1)S (210) 

where t = 4/i 2 — s < 0, u = 0. 

It is a wellknown fact that the t and u channel amplitudes in the case at hand have 
no singularities on the physical sheet of the complex s plane, and, therefore, no imaginary 
parts (see e.g. They are, themselves, imaginary parts of some higher order graphs 

(in fact, of the non-factorizable four-point function of Fig. 14). As a consequence, the 
A4'(t), M.'{u = 0) contributions are small for all t and cannot change the above-discussed 
behaviour. The only point where Ai'if), A4'{u) cause a qualitative change is the elastic 
threshold s = 4/i 2 , t = 0. There the lowest order singular behaviour, equ. (190), that 
was cancelled by the s-channel contribution, equ. (205), is retained, but with a different 
coefficient. We find indeed 



r 



4 

1 loA/C fn\\2 



W 2 (-S ~ 4/i 2 , /i 2 , n 2 ) 



„ 4 l( 2g e g* e E40) + (g$ +^ 2 )E + (0) s2 

W 2 (S ~4/i 2 ,/i 2 ,/i 2 ) V 1- (g g + g*)E+(0) ' 

— > oo for s — > Ajj 2 

where is given by (193) (including the lowest order), and M. '2^2^(0) is given 

by (210) (without lowest order). The E±(0) are just the E± of Section 5 (see (64), (67)). 

Further, whenever the s-channel cross section vanishes (at parity allowed higher pro- 
duction thresholds), its value is changed from zero to a small nonzero number (of order 
(mS) 4 )). 

All the other features of the s-channel scattering cross section remain unchanged. 



XII. SUMMARY 

It is our hope that the discussion of the previous sections has convinced the reader that 
the massive Schwinger model exhibits quite a rich quantum-field theoretic structure and, 
therefore, is worth studying. 

The mass perturbation theory - which was the basic ingredient of our approach - is 
similar to the derivation of a low-energy effective theory in realistic models like QCD. In 
our approach only bosonic states (fermion-antifermion bound states) are present as physical 
states. The (dimensionfull) coupling constant e plays a role similar to Aqcd within our 
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model. In fact, because of the small number of degrees of freedom (Nf = 1) this model 
mimicks, up to a certain extent, the rj' dynamics of QCD. 

A further advantage of the mass perturbation expansion is the fact that - in contrast to 
ordinary perturbation theory - nontrivial phenomena, like instanton vacuum and fermion 
condensates, may be discussed quite straight forwardly (actually, all the computations through- 
out the article are for general vacuum angle 9). 

A first major point in the discussion of the model was the computation of the vacuum 
functional and vacuum energy density. As a consequence, the mass perturbation theory 
could be shown to be IR-finite. 

After deriving the matrix-valued Feynman rules of the mass perturbation theory, we 
could actually give quite a complete description of the physical properties of the model. We 
were able to compute the condensates and susceptibilities and to give an exact description 
of the confinement behaviour. 

Further we succeeded in performing a resummation of the perturbation series for the n- 
point functions with the help of the Dyson-Schwinger equations. These resummed n-point 
functions turned out to be very useful for a further discussion of the physical properties of 
the model. We could infer the whole spectrum of (stable) particles and (unstable) bound 
states from the two-point function. Further we were able to identify all the partial decay 
channels of all bound states. As an illustration, we computed some bound-state masses and 
decay widths. 

At last we discussed scattering processes and found that all the unstable bound states 
turn into resonances of the scattering cross section, as has to be expected. Further, with 
the help of unitarity we were able to identify all the possible final states that may exist for 
the scattering of a given initial state. 

As a result, the following physical picture emerged: there exist two stable particles in 
the theory, namely the Schwinger boson fi and the two-boson bound state Mi- Higher 
(unstable) bound states may be formed out of an arbitrary number of [i and M^. Further, 
these unstable bound states may decay into all combinations of fi and M2 that are possible 
kinematically. (For the special case 9 = the lowest mixed bound state, composed of one 
jj, and one M 2 , is stable, too, because of parity conservation, and is therefore present both 
in final and intermediate states, analogous to our discussion of the M 2 particle; however, 
we treated the generic 9 7^ case throughout the article and, therefore, did not discuss this 
straight-forward generalization in detail.) 

For scattering processes we found that far from all resonances and particle production 
thresholds the scattering cross section is well described by a lowest order computation. 
Whenever the squared momentum s is near a bound-state mass, the scattering cross section 
has a local maximum, i.e. a resonance occurs. Moreover, for all values of s where a new final 
state becomes possible kinematically, the corresponding real particle production threshold 
indeed occurs. 

We want to emphasize that all these features result from our resummed mass perturbation 
theory, and that we did not have to impose further assumptions or use further approxima- 
tions in order to find this physical structure. 

There are (at least) two directions of further study within this field that we believe to be 
important. On one hand, an increase of the number of fermions (Nf > 1) further enriches 
the complexity of the model and makes it possible to mimick the light field dynamics (pions) 
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of QCD, too. There exists some work on the massive multi-flavour Schwinger model that 
uses bosonization and/or semi classical methods ( [^5, 33, 21, 30, 31, 59|). It would be very 
interesting to discuss this model analogously to our discussion in the previous sections; 
however, a direct application of the mass perturbation expansion is not possible there, 
because it fails to be IR-convergent |33|]. Therefore, more elaborate methods should be 
developed for a further study 

On the other hand, the properties we found hold for sufficiently small fermion mass m 
(strong coupling e). Of course, it would be interesting to understand whether and how the 
physical structure of the model changes when the fermion mass is increased. Perhaps a 
discussion similar to this article for the massive Schwinger model within ordinary (electrical 
charge) perturbation theory and a comparison of the two approaches could be a first step 
towards getting more insight into this question. 

In any case, QED 2 remains a fascinating subject of study that will offer further insight 
into general concepts of quantum field theory as well as into some deep problems of QCD. 
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